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A theoretical  and  experimental  study  of  GalnP  and  CdZnSe  visible 
strained  quantum  well  (QW)  lasers  is  presented.  The  strain  in  the  QW 
active  layer  is  either  compressive  or  tensile  and  is  due  to  the  difference  in  lat- 
tice parameter  between  the  QW  and  the  adjacent  sandwiching  barrier  layers. 

Energy  vs.  in-plane  k vector  dispersion  curves  for  strained  QW  active 
layers  are  obtained  using  the  multiband  effective  mass  theory  (or  k • p 
method).  The  carrier  potential  well  structures  are  estimated  using  the  elec- 
tron affinity  rule  for  CdZnSe  and  the  model  solid  theory  for  GalnP  QW  lasers. 
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The  Luttinger-Kohn  Hamiltonian  is  derived  and  the  Schrodinger-like  equa- 
tions  are  solved  nmnerically  by  using  the  finite  difference  method. 

The  spontaneous  emission  polarization  anisotropy  of  OaxIn^.^P 
strained  QW  lasers  for  x=0.3,  0.4  (compressive  strain),  x=0.5  (lattice 
matched),  and  x=0.6  (tensile  strain)  were  measured  and  compared  to  predic- 
tions. The  correlation  between  threshold  current  density  (J^jj)  and  the 
energy  dispersion  curve  of  the  strained  QW  are  also  discussed. 

Incorporating  the  dispersion  curves  calculated  from  the  k • p method 
into  electron-hole  plasma  (EHP)  theory,  the  temperature  dependence  of  of 
CdZnSe  strained  QW  lasers  is  calculated.  Good  agreement  with  experiment 
is  obtained  up  to  300K  provided  a hole-hole  scattering  time  of  60  femtosec- 
onds is  assumed.  Above  300K,  the  measured  J^h  is  much  higher  than  pre- 
dicted, the  discrepancy  being  due  most  likely  to  carrier  leakage. 

To  further  test  the  validity  of  EHP  theory  in  predicting  the  lasing 
behavior  of  CdZnSe  QW  lasers,  dual  wavelength  lasers  (both  CdZnSe  and 
AlGaAs)  were  designed,  fabricated,  and  characterized.  For  the  AlGaAs  QW 
lasers,  dual  wavelength  operation  was  clearly  demonstrated.  For  the 
CdZnSe  QW  lasers,  dual  wavelength  laser  operation  could  not  be  obtained 
although  below  threshold  spectral  data  showed  a double  peak  feature. 
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CHAPTER  I 
INTRODUCTION 


Quantum  confinement  of  carriers  in  a semiconductor  quantum  well 
(QW)  leads  to  the  formation  of  discrete  energy  levels  and  a modified  density 
of  states  function.  A fundamental  advantage  is  that  QW  lasers  deliver  more 
gain  per  injected  carrier  than  conventional  double  heterojimction  (DH)  lasers, 
which  results  in  lower  threshold  currents.  In  addition,  the  spectral  gain 
becomes  quite  anisotropic  with  respect  to  transverse-electric  (TE)  and  trans- 
verse magnetic  (TM)  polarization  when  compared  to  that  of  the  bulk  semicon- 
ductor lasers. 

In  the  early  years  of  semiconductor  lasers,  it  was  thought  that  the  lat- 
tice constants  of  the  active  and  barrier  layer  should  be  matched  in  order  to 
make  a high  quality  heterostructure  diode  laser.  However,  recent  advances 
in  epitaxial  growth  have  allowed  a mismatch  of  the  lattice  constants  between 
the  QW  and  the  barrier  layer.  The  strain  due  to  the  lattice  mismatch  per- 
turbs the  crystal  symmetry  and  leads  to  the  prospect  of  an  artificial  modifica- 
tion of  the  electronic  structure  or  band  structure  of  semiconductor  materials. 
The  combination  of  strain  and  quantiim  confinement  in  the  valence  band  can 
lead  to  substantially  more  favorable  energy  dispersion  curves  for  laser  action 
than  those  existing  in  the  imstrained  semiconductor  crystal. 
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In  recent  years,  the  conjecture  [Yabl  86,  Adam  86]  that  strained  active 
layers  in  QW  lasers  should  lead  to  improved  performance  compared  with  lat- 
tice matched  lasers  has  been  verified  with  visible  lasers  by  experiment.  In 
fact,  a strained  QW  design  was  utilized  in  fabricating  the  world’s  first  blue- 
green  diode  lasers  [Haas  91].  Moreover,  it  has  been  shown  that  strained  QW 
lasers  can  have  reliability  as  well  as  performance  characteristics  that  are 
substantially  superior  to  unstrained  devices.  In  this  study,  GalnP  (red)  and 
CdZnSe  (blue-green)  strained  QW  lasers  are  analyzed.  The  applications  and 
the  fimdamental  characteristics  to  be  analyzed  are  listed  in  Table  1.1. 

In  semiconductor  lasers,  the  condition  required  for  positive  gain  is 

Eg  < hv  < Efj,  - Efy  (1.1) 

where  Eg  is  the  band  gap  energy,  hv  is  the  photon  energy,  and  Efj.(Efy)  is  the 
quasi-Fermi  level  for  the  conduction  (valence)  band.  The  condition  in  Eq. 
(1.1)  was  first  pointed  out  by  Bernard  and  Duraffourg  [Bern  61].  When  the 
quasi-Fermi  level  separation  is  equal  to  the  band  gap  energy,  the  material 
will  become  transparent  for  the  photon  energy  equal  to  the  band  gap.  The 
electron  and  hole  density  that  is  required  to  provide  this  separation  is  known 
as  the  transparency  carrier  density,  N^j.,  and  its  magnitude  is  fimdamentally 
related  to  the  densities  of  states  of  conduction  ( p ) and  valence  band  ( p ) of  a 
given  material.  The  transparency  condition  is  illustrated  for  my=m^  and 
mv>m(.  where  mc(my)  is  the  effective  mass  of  the  conduction  (valence)  band  in 
Fig.  1.1.  In  the  latter  case,  the  asymmetry  of  the  density  of  state  shifts  both 
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Table  1.1 

The  Characteristics  of  the  Red  eind  the  Blue-Green  Lasers 


Red  Lasers 

Blue-Green  Lasers 

Material 

GaxIni-xP/ 

Cd.2Zn.8Se/ 

(QW/Barrier) 

(Alo.6Gao.4)o.5lno.5P 

ZnS.oeSe.94 

Wavelength 

707nm  (x=0.3) 
630nm  (x=0.6) 

496nm  (85K) 
516nin  (294K) 

Applications 

• Optical  bar  code 
scanning 

• Optical  Recording 

• He-Ne  laser 
substitutes 

• Laser  printers 

• Optical  memory 

• Compact  disk  player 
optical  source 

• Under  sea 
communications 

• Color  display 

Characteristics 

s 

Polarization  Anisotropy 

Temperature 
Dependence  of  Jth 
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Ef, 


. E 


fc 


Efy 


Figure  1.1  Band  edge  tranparency  condition  (Efg-Efy=Eg)  for  (a)  m(.=my 

and  (b) 
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Fermi  fimctions  toward  the  conduction  band  in  order  to  maintain  equal  num- 
bers of  carriers  in  both  bands.  The  result  is  a larger  N^j.  than  in  the  symmet- 
ric case. 

In  most  semiconductor  lasers,  there  is  a serious  asymmetry  between 
the  light  conduction  band  effective  mass  and  the  very  heavy  valence  band 
effective  mass.  Ideally  both  masses  should  be  as  light  as  possible.  The  den- 
sity of  states  woiold  be  very  small  and  the  carrier  density  required  for  the 
transparency  would  be  minimized.  Due  to  the  heavy  valence  band  effective 
mass,  the  hole  quasi-Fermi  level  is  above  the  top  of  the  valence  band  and  the 
hole  distribution  is  classical.  This  situation  is  illustrated  in  Fig.  1.2(a)  and 
the  more  ideal  situation  is  illustrated  in  Fig.  1.2(b).  Due  to  the  effective 
mass  asymmetry  of  the  typical  III-V  semiconductor  as  illustrated  in  Fig. 
1.2(a),  the  quasi-Fermi  levels  are  offset  from  the  symmetrical  position  by  an 
energy  5.  This  requires  an  excess  population  of  carriers  as  is  evident  by  com- 
paring the  electron  Fermi  level  in  Fig.  1.2(a)  and  (b).  When  both  hole  and 
electron  masses  are  light  as  in  Fig.  1.2(b),  the  transparency  condition  can  be 
achieved  at  lower  carrier  injection. 

The  valence  band  effective  mass  can  be  reduced  significantly  by 
applying  strain  in  the  QW  layer.  The  strain  splits  the  degeneracy  of  the 
heavy  and  light  hole  states,  accessing  a wide  range  of  subband  structures, 
including  the  possibility  of  the  highest  valence  subband  being  light  hole-like, 
of  significant  benefit  for  semiconductor  lasers  in  terms  of  the  threshold 
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Figure  1.2  The  E vs.  k band  structure  diagram  of  (a)  a typical  III-V  semi- 
conductor and  (b)  an  idealized  strained  QW. 
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current  density.  Moreover,  the  Auger  recombination  and  the  intervalence 
band  absorption,  which  are  the  major  loss  mechanisms  in  long  wavelength 
lasers,  can  be  reduced  [Yabl  86,  Wang  93]. 

Another  change  accompanied  by  the  strain  effect  involves  the  polariza- 
tion anisotropy.  The  transition  probability  between  the  conduction  band  and 

« 

heavy  hole  band  is  greatest  for  TE  polarization.  Conversely,  TM  polarization 
is  strongest  for  the  conduction  band  and  the  light  hole  band  transition.  For 
most  unstrained  semiconductor  lasers,  the  heavy  hole  band  is  the  ground 
state.  Therefore,  the  heavy  hole  population  is  most  easily  inverted  and  the 
gain  polarization  asymmetry  of  the  heavy  hole  transition  favors  the  TE  mode 
over  TM  mode.  If  the  quantiun  well  active  layer  is  strained,  the  polarization 
anisotropy  is  modified.  For  instance,  Bour  et  al  [Bour  93]  fabricated  and 
analyzed  red  light  emitting  GaxIni.xP/  (Alo,eGao.4)o.5lno.5P  strained  QW 
lasers  grown  on  GaAs  substrates.  For  the  Gallium  (Ga)  mole  fraction  (x)  less 
than  0.5,  the  quantum  well  layer  is  compressively  strained  and  the  heavy 
hole  is  the  ground  state.  In  this  case,  TE  mode  emission  dominates.  On  the 
other  hand,  for  the  tensile  strain  (x  greater  than  0.5),  the  light  hole  is  the 
ground  state  and  the  emission  is  predominantly  TM.  In  either  case,  the 
polarization  anisotropy  is  important  in  order  to  predict  the  diode  laser  perfor- 
mance. 

In  this  study,  a theoretical  and  experimental  study  of  GaxIn^.xP/ 
(Alo,eOao.4)o.5lno  5P  (red)  and  Cdo  2Zno.8Se/ZnSo.6Se  (blue-green)  strained 
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QW  lasers  is  presented.  For  OaxIn^.^P  QW  lasers,  polarization  anisotropy  of 
the  compressive  as  well  as  the  tensile  strain  is  analyzed.  The  electron  hole 
plasma  (EHP)  model,  that  optical  gain  is  derived  from  k-conserved  stimu- 
lated processes  between  free  electrons  and  holes  in  quantum  well  subbands, 
is  used  for  CdZnSe  QW  lasers  to  predict  the  temperature  dependence  of  the 
threshold  current  density  and  compared  to  the  recently  demonstrated  [Gain 
93]  experimental  results. 

In  chapter  II,  the  energy  vs.  in-plane  dispersion  curves  for  red  and 
blue-green  lasers  are  calculated  using  the  multiband  effective  mass  theory. 
The  Luttinger-Kohn  Hamiltonians  are  derived  in  a simplified  form  and  the 
Schrodinger  equations  are  solved  using  the  finite  difference  method  (FDM). 
The  theories  for  a heterojimction  band  lineup  are  introduced  and  the  strain 
effects  on  the  band  offsets  are  discussed.  The  carrier  potential  wells  for 
GaxIni_xP/(Alo  6Gao.4)o.5lno.5P  and  Cdo  ^ZnQ  sSe/ZnSo  gSe  strained  QW  lasers 
are  derived  with  the  model  solid  theory  and  the  electron  affinity  rule,  respec- 
tively. 

In  chapter  III,  the  polarization  anisotropy  of  GaxIn^.xP/ 
(Alo.6Gao.4)o.5lno.5P  strained  QW  lasers  for  x=0.3,  0.4  (compressive  strain), 
x=0.5  (lattice  matched),  and  x=0.6  (tensile  strain)  is  described  in  detail. 
Based  on  the  energy  dispersion  curves  derived  in  chapter  II,  the  peak  emis- 
sion differences  between  the  TE  £md  TM  polarization  are  calculated  and  com- 
pared to  the  experiment.  The  threshold  current  densities  are  measured  and 
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discussed  in  terms  of  the  in-plane  effective  masses  and  quasi-Fermi  levels. 

In  chapter  IV,  the  temperature  dependence  of  threshold  current  den- 
sity of  Cdo.2Zno,8Se/ZnSo.6Se  is  predicted  and  compared  to  the  experiment. 
A modified  Lorentzian  line  shape  fimction  is  used  to  include  the  intraband 
relaxation  process  caused  by  carrier  scattering.  The  spontaneous  emission 
rate  and  the  spectral  gain  are  calculated  based  on  the  energy  dispersion 
curve  derived  in  chapter  II.  The  peak  gain  (g)  vs.  the  current  density  (J) 
curves  at  various  temperatures  are  derived  by  relating  g and  J at  each  carrier 
density.  The  temperature  and  the  carrier  density  are  specified  first,  and 
then  the  spontaneous  emission  rate  and  the  spectral  gain  are  calculated. 
The  corresponding  J is  calculated  by  integrating  the  spontaneous  emission 
rate  and  g is  determined  from  the  maximum  of  the  spectral  gain.  The 
threshold  current  densities  at  various  temperatures  are  determined  from  the 
g-J  curve  by  estimating  the  threshold  gain. 

There  have  been  reports  [Ding  90,  Ding  92]  that  the  optical  gain  in 
CdZnSe  QW  lasers  is  not  derived  from  stimulated  processes  between  free 
electrons  and  holes  but  between  bound  electron  and  holes  or  excitons. 
However,  the  predicted  temperature  dependence  of  threshold  current  density 
of  Cdo.2Zno  gSe/ZnSo.eSe  QW  lasers  is  in  good  agreement  with  experiment 
without  any  excitonic  effect  consideration.  To  further  test  the  validity  of 
EHP  theory  in  predicting  the  lasing  behavior  of  CdZnSe  QW  lasers,  dual 
wavelength  lasers  (both  CdZnSe  and  AlGaAs)  were  designed,  fabricated,  and 
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characterized  as  discussed  in  chapter  V.  For  the  AlGaAs  QW  lasers,  dual 

wavelength  operation  was  clearly  demonstrated.  For  the  CdZnSe  QW  lasers, 

the  spectral  data  appeared  to  show  dual  wavelength  operation  although  with 

« 

considerably  reduced  resolution. 


CHAPTER  II 

THE  MULTIBAND  EFFECTIVE  MASS  THEORY  AND 

ENERGY  DISPERSION 


2.1  Introduction 

The  multiband  effective  mass  theory,  often  called  k • p method,  has 
been  widely  used  to  treat  problems  involving  valence  band  and  strongly  non- 
parabolic bands  in  narrow  gap  semiconductors  [Dres  55,  Lutt  55,  Kane  57, 
Piku  60].  In  recent  times,  it  has  been  used  to  describe  quantum  wells  and 
superlattices  [Whit  81,  Bast  81,  Alta  85]  and  regarded  as  the  single  most 
important  theoretical  procedure  for  predicting  and  analyzing  the  energy  band 
structures  of  semiconductors. 

In  order  to  calculate  the  energy  vs.  in-plane  k vector  dispersion  of 
GaxIni.xP/(Alo.6Gao,4)o.5lno.5P  and  Cd0.2Zn0.8Se/  ZnS0  6 QW  lasers  with  the 
multiband  effective  mass  theory,  one  must  first  define  the  carrier 
confinement,  or  potential  well  structure  of  these  material  systems.  The 
carrier  confinement  can  be  determined  either  by  experiment  such  as 
magneto-optical  spectroscopy  or  by  theory  such  as  the  model  solid  theory 
[Van  86,  Van  88,  Van  89]  or  the  electron  affinity  rule.  Since  the  experimental 
data  are  not  known  for  these  materials,  the  carrier  confinement  was 
estimated  theoretically.  Moreover,  the  lattice  constants  of  the  QW  layers  of 
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these  structure  are  different  from  those  of  the  barrier  layers,  and  hence,  the 
QW  layers  are  strained.  Since  the  strain  perturbs  the  band  structure,  it 


In  section  2.2,  the  band  effective  mass  theory  is  developed  from  the 
fundamental  Bloch  theorem  and  in  section  2.3,  the  carrier  potential  wells  for 

GaxIni.xP/(Alo.6Gao.4)o.5lno.5P  and  Cdo.2Zno.8Se/  ZnSo.o  QW  laser  structures 

are  derived  with  the  strain  effect  considered.  The  energy  vs.  in-plane  k vec- 

« 

tor  dispersion  curves  are  then  calculated  using  the  multiband  effective  mass 
theory  in  section  2.4. 

2.2  The  Multiband  Effective  Mass  Theory  and  Valence  Band  Mixing 

The  Bloch  theorem  states  that  the  eigenfrmctions  of  the  wave  equation 
for  a periodic  potential  V(r)  are  of  the  form  of  the  product  of  a plane  wave 

ik  * r , 

e times  a frmction  u . (r)  with  the  periodicity  of  the  crystal.  Hence,  the 
eigenvalue  equation  for  the  energy  E in  band  n is . 


shotdd  be  taken  into  account  in  the  carrier  confinement  calculation  section. 


(2.1) 


(2.2) 


where  u . (r)  is  the  normalized  Bloch  function  with  the  periodicity  of  V(r) 
and  wave  vector  k.  Assume  that  the  Bloch  function  is  a slow  varying 
frmction  of  k in  a localized  state  such  as  in  quantum  wells.  Then  the  Bloch 
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ftmction  with  arbitrary  k is  approximately  where  UjjQ(r)  is 

the  Bloch  ftmction  at  the  band  edge  k=0.  Since  t)/  . ( r)  is  complete  for 
functions  having  the  periodicity  of  V(r)  for  any  given  k [Lutt  55],  a localized 
wave  function  (r)  may  be  expressed  in  terms  of  the  linear  combination  of 
\|/  k (r)  using  an  arbitrary  set  of  coefficients  A^Ck). 


'I'nk  C)  = SA„0‘)e“‘-\,(r) 

k 

k 


(2.3) 


where  F„(r)  is  £m  envelope  function.  The  description  used  in  Eq.  (2.3)  is 
known  as  the  envelope  function  approximation.  By  combining  Eq.  (2.1) 
through  Eq.  (2.3),  the  matrix  element  of  the  Hamiltonian  H can  be  derived  as 
follows: 


(n’,k'|H|n,k>  = 


2 9^ 


^nn’^k,  k’  m 


k • p .5,  , . 

k,  k 


(2.4) 


o 


<\,olPlVo> 


= -i  R (u 


(2.5) 


The  detailed  derivation  of  Eq.  (2.4)  and  Eq.  (2.5)  is  shown  in  the 
Appendix.  The  last  term  in  Eq.  (2.4),  as  well  as  the  usually  very  small  term 
in  k2  , can  be  treated  as  a perturbation  of  the  normal  crystal  Hamiltonian. 

The  Hamiltonian  in  matrix  form  can  be  constructed  using  Eq.  (2.4)  and 
Eq.  (2.5).  For  simplicity,  it  is  assumed  that  the  conduction  and  valence 
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bands  are  decoupled,  which  is  in  fact  a good  approximation  in  a wide  band 
gap  material.  The  matrix  is  then  derived  in  terms  of  |J,  mj)  representation, 
for  J=3/2  by  considering  a heavy  hole  and  a light  hole  band  with  spin  up  and 
down  in  each  state  and  shown  in  Table  2.1.  The  four  by  four  matrix  can  be 
greatly  simplified  using  the  method  first  suggested  by  Kane  [Kane  57]  and 
later  refined  by  Broido  et  al  [Broi  85,  Broi  86],  Altarelli  et  al  [Alta  85],  and 
Ahn  et  al  [Ahn  88].  They  pointed  out  that  a imitary  transformation  with  a 
proper  choice  of  basis  set  can  decouple  the  four  by  four  matrix  into  a pair  of 
identical  two  by  two  matrix  systems. 

With  this  simplification,  the  Schrodinger  equation  for  heavy  and  light 
holes  is  expressed  as 


H,W 

Ph 

EhPh 

Wt  Hj 

Fl 

EjFi 

(2.6) 


where  Hjj,  Hj  are  the  Hamiltonian  of  the  heavy  and  light  holes,  W and  W”  are 
the  coupling  terms  and  Fh  and  Fj  are  the  envelope  fimctions  of  the  heavy  and 
the  light  holes.  Assuming  the  plane  wave  solutions,  Hjj  and  can  be  writ- 
ten as 


Hj  = (Yi  + 2y^)  kl  + (Yj  - Y2)  kj  + Vj 


(2.7) 


where  Yj's  the  Luttinger  parameters  and  Vj  are  the  potential  well 
depth  of  the  heavy  and  fight  holes.  In  the  above  equation,  and  kj;  are  the  k 
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Table  2.1 

The  Luttinger-Kohn  Matrix  for  Heavy  and  Light  Hole  Bands 


1- 

'2’  2' 

i--i> 

7’  2' 

-> 
T 2' 

|3  J) 
'r  1' 

3 3 
2’ 2 

P+Q 

R 

-S 

0 

|3  J) 

7’  2' 

R* 

P-Q 

0 

s 

1^  -> 
'2’  2' 

-S* 

0 

P-Q 

R 

|3  J) 
'2’  2' 

0 

s* 

R* 

P+Q 

P = r 


1 

2 


Q = ^ 


1 

2 


y.fk^  + k^  + k^l 

* IV  X y zj 


2 .2  ,2 
+ k 
y z 


Yjl  k + k + k 


s = S 


v“>o/ 


( 

R = - 

h 

Vak.  (kx  - iky) 
X (V2 + y,)  (k 


•1  73. 


Y2)  (k^  + ik  ) 
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vector  components  with  directed  along  the  epi  growth  direction  and  kj. 
directed  along  the  epi  layer.  The  coupling  terms  W and  W"*"  for  the  {100} 
plane  are  defined  as  [Zory  93] 


W = 73kj(Y2kj-i2Y3k^) 
Wt  = V3k^(Y2kj  + i2Y3k^)  . 


For  kj.=0,  the  coupling  term  vanishes  so  we  can  solve  the  coupled 
Schrodinger  like  equations  for  the  valence  band  separately,  but  for  finite  kj., 
we  need  to  solve  them  simultaneously.  It  is  noted  that  from  Eq.  (2.7),  the 
Luttinger  parameter  and  the  effective  masses  can  be  related  for  k^=0  as  fol- 
lows [Zory  93]: 


Yi-2Yj 


7i+2y2  2m, I, 


(2.9) 


2.3  Carrier  Potential  Well 

In  order  to  calculate  the  energy  dispersion  of  given  QW  structures,  it  is 
essential  to  know  the  carrier  potential  well  structure,  or  potential  well  depth. 
In  this  section,  the  quantum  well  structures  of  GaxIni.xP/CAlo.eGao  4)o.5lno,5P 
and  Cdo.2Zno,8Se/ZnSo,o6^®  estimated.  In  general,  determining  the  band 
lineups  in  a heterostructure  is  a very  complicated  problem.  So  it  is  hard  to 
predict  the  band  offsets  with  great  accuracy  even  with  the  best  theory  known. 
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In  other  words,  the  exact  band  offset  can  be  determined  only  using  experi- 
mental techniques  such  as  magneto-optical  spectroscopy.  However,  imfortu- 
nately,  the  experimental  data  for  the  band  offsets  of  these  material  systems 
are  not  known  yet.  In  this  section,  two  methods  for  estimating  the  carrier 
confinement  will  be  introduced:  the  electron  affinity  rule  which  is  the  simplest 
model  and  the  model  solid  theory  by  Van  de  Walle  and  Martin  [VAN  86] 
which  is  nowadays  widely  accepted  and  supported  by  experiment. 


2.3.1  The  Model  Solid  Theory 

Recently,  Van  de  Walle  et  al  [Van  88]  have  reported  a series  of  calcula- 
tions for  a number  of  heterojunction  interfaces.  In  that  work,  the  large  num- 
ber of  interfaces  is  treated  in  a uniform  way.  In  the  model  solid  theory,  the 
common  reference  energy  level  is  accomplished  by  modeling  the  solid  as  a 
superposition  of  neutral  atoms.  In  each  atom,  the  electrostatic  potential  is 
rigorously  defined  with  respect  to  the  vacuum  level.  Therefore,  the  average 
electrostatic  potential  (Egy)  is  specified  on  the  absolute  energy  scale  by  super- 
position. The  relative  positions  of  valence  band  (Ey)  and  conduction  band 

«• 

(Eg)  edge  with  respect  to  the  average  potential  are  simply  given  by 


V,  av 


+ 


(2.10) 


where  Aq  is  the  spin-orbit  splitting  and  Eg  is  the  energy  band  gap.  The 
average  potential  itself  contains  no  information.  But  when  two  different 
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semiconductors  consist  of  a heterojimction,  we  can  define  the  band  offset  by 
comparing  the  average  potential  (and  hence,  E^'s  and  Ey's)  of  each  material. 
The  band  offset  parameters  of  GaxInx.xP/(Alo.6Grao.4)o.5lno.5P  for  x=0.3,  0.4, 
0.5,  and  0.6  are  summarized  in  Table  2.2. 

2.3.2  The  Electron  Affinity  Rule 

Although  the  model  solid  theory  is  regarded  as  significantly  more  accu- 
rate than  the  earlier  ones  for  heterostructure  band  lineups,  the  average  elec- 
trostatic potential  (E^y)  of  CdSe  (which  is  needed  for  interpolation  for 
Cdo.2Zno.8Se)  is  not  known  yet  because  of  the  lack  of  material  parameters  for 
cubic  CdSe  (most  hexagonal  CdSe  parameters  are  known).  Hence,  the  sim- 
ple electron  affinity  rule  was  used  for  the  CdZnSe  system. 

To  define  the  band  lineups  of  the  Cdo.2Zno  8Se/ZnSo.eSeo.94  heterojunc- 
tion, the  conduction  and  valence  band  offset  in  the  absence  of  strain  was  esti- 
mated using  the  electron  affinity  rule.  The  vacuum  level  was  used  for  a 
reference  and  the  relative  positions  of  the  conduction  band  edges  of 
Cdo.2Zno.8Se  and  ZnSo.eSeo.94  were  determined  from  the  electron  affinity  of 
each  material,  that  is,  the  energy  to  remove  an  electron  from  the  conduction 
minimmn  to  the  vacuum  level.  Then  the  valence  band  offset  was  derived  by 
subtracting  the  energy  band  gap  from  the  conduction  band  edge  of  each  mate- 
rial. The  electron  affinity  is  4.079eV  for  ZnSo.6Seo.94  and  4.262eV  for 
Cdo.2Zno  8Se  as  interpolated  from  the  host  materials  (ZnS,  ZnSe  and  CdSe). 
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Table  2.2 

GaxInj.xP  Parameters 


Ga^Ini.xP 

®v,av 

Ao 

Eg 

Ev 

Ec 

x=0  (InP) 

-7.04 

0.11 

1.349 

-7.003 

-5.654 

x=0.3 

-7.15 

0.10 

1.624 

-7.117 

-5.493 

o 

II 

X 

-7.19 

0.10 

1.746 

-7.157 

-5.411 

• 

o 

II 

X 

-7.23 

0.10 

1.883 

-7.197 

-5.314 

x=0.6 

-7.26 

0.108 

2.035 

-7.224 

-5.189 

x=l  (GaP) 

-7.40 

0.08 

2.797 

-7.373 

-4.576 

(Al  5Ga_4)  jln^P 

-7.44 

0.13 

2.31 

-7.397 

-5.087 
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The  problem  of  this  approach  is  that  the  electron  affinity  is  not  neces- 
sarily the  property  of  the  heterojunction  interface  [Capa  87].  Another  prob- 
lem with  the  electron  affinity  rule  is  that  there  may  be  charge  transfer  at  the 
interface,  which  would  alter  the  band  lineup.  So  it  should  be  noted  that  our 
results  will  be  just  a rough  approximation. 


2.3.3  The  Strain  Effect 

When  the  lattice  constant  of  QW  material  is  different  from  the  barrier 
material,  the  QW  layer  is  subject  to  either  compressive  or  tensile  strain  as 
illustrated  in  Fig.  2.3.  The  strain  reduces  the  crystal  symmetry  and  results 
in  the  energy  shifts  of  the  conduction  and  valence  band.  As  a consequence, 
the  energy  vs.  in-plane  k vector  dispersion  varies  and  hence  the  in-plane 

9 

effective  masses  are  changed.  The  energy  shifts  due  to  the  hydrostatic  strain 
are  expressed  as  [Poll  68,  Gavi  70,  Chan  77] 


AE  (hydrostatic)  = 

c ^ L2 

AEy  (hydrostatic)  = 


for  conduction  band 


for  valence  band 


(2.11) 


where  ac(a„)  are  hydrostatic  deformation  potential  for  the  conduction 


(valence)  band,  and 


AQ  . 

H 


is  the  fractional  volume  change.  For  (001)  sub- 


strates, AQ/Q  = £v-ir  ^vv  ^11  where  £••  are  the  components  of  the  strain 

A A yy  Z.Z.  Ij 

tensor  defined  ase  =e  =(a./a  )-l  for  the  direction  parallel  to  the 

XX  yy  CD  qw^  ^ 

surface,  where  at  (aqw)  is  the  lattice  constant  of  the  cladding  (quantum  well) 

2Ci2 


layer  and  = 


C„  XX 


e for  the  direction  perpendicular  to  the  surface 
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Barrier 


Barrier 


Figure  2.3  Schematics  of  strain  effect  for  (a)  compressive  and  (b)  tensile 


strain. 
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where  Cy  are  the  elastic  moduli.  The  shear  strain  splits  the  degeneracy  of 
the  valence  band  at  F peak.  For  the  strain  along  the  (001)  plane,  the  energy 
shifts  of  heavy  and  light  hole  states  are  expressed  as  [Poll  68,  Gavi  70,  Chan 
77] 


AE^  (shear)  = 
AEj(shear)  = -^  - 


1 (SEooi)^ 


(2.12) 


with  SEqqj  = 2b  (e  - e ) where  b is  the  shear  deformation  potential.  The 
parameters  for  band  lineups  are  summarized  in  Table  2.3. 


2.3.4  The  Potential  Well  Structure  of  GalnP  and  CdZnSe  QW  Lasers 

Based  on  the  methods  discussed  in  the  previous  sections,  the  potential 
well  structures  of  GaxIn^.xPAAlo  6Gao.4)o.5lno.5P  and  Cdo  ^ZnQ  gSe/ZnSo.oeSe 
quantum  well,  shown  in  Fig.  2.4  and  Fig.  2.5,  are  calculated  in  two  steps; 
First,  the  band  offset  is  estimated  in  the  absence  of  the  strain  effect  with  the 
model  solid  theory  for  Ga^In^.^PAAlo  0Gao,4)o.5lno.5P  quantiim  wells  and  with 
electron  affinity  for  Cd0.2Zn0.8Se/ZnS0.06Se  quantum  wells.  Then  the  strain 
effect  is  taken  into  account  for  each  quantum  well  and  finally  the  carrier  con- 
finements are  derived.  The  results  are  shown  in  Fig.  2.6,  Fig.  2.7,  and  Fig.  2.8 
schematically  for  GaxIni_xP/(Alo.6Gao.4)o.5lno.5P  for  x=0.3(compressive  strain) 
and  x=0.6  (tensile  strain)  and  Cd0.2Zn0.8Se/ZnS0.06Se  quantum  well  lasers. 

For  GaxIn4.xP/(Alo.6Gao.4)o.5lno.5P  quantum  wells,  an  important 
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Table  2.3 

GaxInj.xP  Material  Parameters 


InP 

x=0 

CO 

• 

o 

II 

X 

X 

II 

o 

x=0.5 

X 

II 

o 

GaP 

X=1 

ao 

5.8686 

5.7434 

5.7016 

5.6599 

5.6182 

5.4512 

Cll 

10.22 

11.46 

11.86 

12.25 

12.66 

14.12 

C12 

5.76 

5.92 

5.97 

6.02 

6.07 

6.25 

ay 

1.27 

1.41 

1.45 

1.49  - 

1.53 

1.70 

ac 

-5.04 

-5.72 

-5.92 

-6.12 

-6.33 

-7.14 

a 

-6.31 

-7.12 

-7.37 

-7.62 

-7.87 

-8.84 

b 

-2.0 

-1.94 

-1.92 

-1.90 

-1.88 

-1.80 

™hh 

0.45 

0.464 

0.462 

0.460 

0.458 

0.47 

mih 

0.12 

0.135 

0.14 

0.145 

0.15 

0.17 
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Figure  2.4  The  epitaxial  layer  structure  of  Ga^Iui.^P  QW  lasers. 
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Figure  2.5 


The  epitaxial  layer  structure  of  CdZnSe  QW  lasers. 
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Figure  2.6  Band  line-ups  of  Gao.sino  7P/(Alo  6Gao.4)o.5lno.5P  QW  struc- 
ture; (a)  bulk  (b)  strain  effect  (c)  carrier  confinement  with 
strain. 
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Figure  2.7  Band  line-ups  of  Gao  6lno.4P/(Alo.6Gao.4)o.5lno.5P  QW  struc- 
ture; (a)  bulk  (b)  strain  effect  (c)  carrier  confinement  with 
strain. 
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Figure  2.8  Band  line-ups  of  Cdo.2Zno.8Se/ZnSo.oeSe  QW  structure; 

(a)  bulk  (b)  strain  effect  (c)  carrier  confinement  with  strain. 
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feature  can  be  seen.  The  heavy  hole  is  the  ground  state  for  compressive 
strain  (x=0,3)  while  the  light  hole  is  the  grormd  state  for  tensile  strain 
(x=0,6).  These  differences  lead  to  polarization  anisotropy  and  will  be 
discussed  in  detail  in  chapter  III.  For  CdZnSe  quantiun  wells,  the  potential 
well  depth  of  the  heavy  hole  state  is  72meV,  but  the  light  hole  state  is  almost 
lined  up  with  the  barrier  layer.  So  it  is  expected  that  most  transitions  occur 
between  conduction  electrons  and  relatively  well  confined  heavy  hole  states. 


2.4  Energy  Dispersion 

Based  on  the  band  diagram  under  strain  effect  derived  in  the  previous 
section,  we  now  derive  the  energy  vs.  in-plane  k vector  dispersion.  We 
assume  that  the  conduction  band  is  parabolic  and  take  the  band  mixing  effect 
into  account  for  the  valence  band.  The  multiband  effective  mass  theory  was 
used  to  calctJate  the  band  structure  and  electron  (hole)  wave  fimctions  of  a 
quantum  well. 

For  the  conduction  band,  the  subband  structure  of  a quantmn  well  is 
obtained  by  solving  the  Schrodinger-like  equation  using  the  envelope  function 
approximation. 

H F,  = E (2.13) 

c c c c 


where  is  the  hamiltonian  for  the  conduction  band  simply  given  by 


c 2m  dz^  ^ 


(2.14) 


where  F^  is  the  envelope  function  of  the  conduction  band,  m^.  is  the 
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conduction  band  effective  mass,  and  Vj.  is  the  potential  well  depth  for  the 
conduction  band  electron.  The  assumption  that  the  interaction  between 
energy  subbands  is  weak  makes  the  calculation  of  the  conduction  band 
electron  structure  simple.  Since  this  does  not  hold  in  the  valence  band,  the 
k • p method  was  used  for  the  valence  band  and  band  mixing  was  taken  into 
account. 

The  finite  difference  method  (FDM)  was  used  for  the  numerical  calcu- 
lation of  the  energy  dispersion.  This  method  is  particiilarly  well  smted  to 
simple  device  geometries  and  has  been  used  extensively  to  model  one  and  two 
dimensional  rectangular  geometry  (planar)  devices  such  as  quantum  well 
lasers.  The  FDM  produces  solutions  for  the  physical  variables  as  discretized 
values  at  specific  nodes  contained  within  a mesh  superimposed  on  the  solu- 
tion domain.  The  finite  difference  equations  are  derived  from  truncated  Tay- 
lor series  for  continuous  function  f(x+Ax)  and  fix- Ax)  [Lim  92].  For 
simplicity,  the  uniform  meshes  were  used  in  the  calculation.  The  quantum 
well  region  was  divided  into  20  to  40  regions,  depending  on  the  thickness, 
and  about  200  mesh  points  in  the  barrier  layer  were  included  for  a better  con- 
vergence.  Then  the  Schrodinger  equations  Eq.  (2.6)  and  Eq.  (2.13)  were  dis- 
cretized at  each  node  as  schematically  illustrated  in  Fig.  2.9.  The  energy 
dispersion  for  GaxIni.xP/(Alo,eGao.4)o.5lno.5P  and  Cdo  ^Zng  gSe/  ZnSo  e QW 
lasers  are  shown  in  Fig.  2.10,  Fig.  2.11,  andFig.  2.12. 
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Figure  2.9  Mesh  points  and  discretized  function  of  a QW  structure 
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Figure  2,10  The  energy  vs.  in-plane  k vector  dispersion  of  Oaxln^.x?/ 

(Alo.6Gao.4)o.5lno.5P.  (a)  x=0.3  and  (b)  x=0,4. 
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Figure  2.11  The  energy  vs.  in-plane  k vector  dispersion  of  GaxInj.xP/ 

(Alo.6Gao.4)o.5lno.5P.  (a)  x=0.5  and  (b)  x=0.6. 
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Figure  2.12  The  energy  vs,  in-plane  k vector  dispersion  of  Cdo.2Zno.8Se/ 

ZnSo.6  QW  lasers. 
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2.5  Summary 

The  energy  vs.  in-plane  k vector  dispersions  of  GalnP  and  CdZnSe 
quantum  wells  were  derived.  For  the  conduction  hand,  the  parabolic  approx- 
imation was  used  neglecting  the  interaction  between  the  energy  subbands. 
For  the  valence  band,  the  multiband  effective  mass  theory  was  used  and  the 
band  mixing  effect  was  taken  into  accotmt.  In  order  to  define  the  carrier 
potential  well  structures,  the  electron  affinity  rule  and  the  model  solid  theory 
were  used.  The  strain  on  the  QW  layer  perturbed  the  crystal  symmetry  and 
resulted  in  energy  shifts  of  the  conduction  emd  valence  bands.  The  streiin  is 
either  compressive  or  tensile  and  the  energy  gap  is  increased  for  a compres- 
sive strain  and  reduced  for  a tensile  strain  by  a hydrostatic  deformation. 
The  shear  strain  splits  degeneracy  of  the  valence  band  at  the  band  edge. 
With  the  strain  effect  taken  into  accoimt,  it  was  shown  that  the  valence  band 
ground  state  of  OaxIn^.^P/fAlo  eGao  4)o.5lno,5P  QW  lasers  can  be  either  a 
heavy  hole  state  or  light  hole  state  depending  on  the  GaP  alloy  content.  For 
CdZnSe  QW  lasers  whose  active  layer  is  compressively  strained,  the  light 
hole  state  was  almost  lined  up  with  the  barrier  layer  while  the  heavy  hole 
state  was  well  confined.  The  energy  dispersion  curves  for  GalnP  and 
CdZnSe  QW  lasers  were  numerically  calciilated  using  the  finite  difference 


method. 


CHAPTER  III 

POLARIZATION  ANISOTROPY  OF  RED  LASERS 

3.1  Introduction 

Red  diode  lasers  with  peak  wavelength  around  680  nm  are  of  interest 
for  numerous  application  areas.  Important  uses  include  laser  pointers,  solid 
state  laser  pumping,  medicine,  printing,  bar  code  scanning,  and  optical  data 
storage.  It  is  known  that  such  lasers  are  most  successfully  made  with  GalnP 
quantum  well  lasers  with  AlGalnP  barrier  layer.  However,  depends  on  the 
GaP  alloy  content  in  the  active  layer,  the  quantum  well  is  subject  to  the 
strain. 

Shear  strain  perturbs  the  valence  band  structure  of  quantum  wells, 
and  are  important  in  determining  the  threshold  current  density,  wavelength, 
and  polarization  state  of  quantum  well  diode  laser.  In  fact,  quantum  well 
lasers  most  often  oscillate  in  the  transverse  electric(TE)  polarization,  because 
the  heavy  hole(hh)  band  is  the  groimd  state  for  both  lattice  matched  and  com- 
pressive strained  quantum  wells  [Chin  88,  Zory  93].  Therefore,  the  hh  popu- 
lation is  most  easily  inverted,  and  the  polarization  anisotropy  of  the  hh 
transition’s  gain  favors  the  TE  mode.  Conversely,  transverse  magnetic(TM) 
mode  oscillation  prevails  in  some  tensile  strained  quantum  well  lasers,  where 
the  light  hole  (Ih)  band  becomes  the  ground  state  through  a strain  induced 
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reversal  of  the  Ih  £ind  hh  band  edges  [Welc  91,  Boer  90,  Thij  91].  In  this 
chapter,  the  polarization  dependent  effects  and  the  in-plane  spontaneous 
emission  of  80A  GaxIni.xP/  (Alo.6Gao.4)o.5lno.5P  quantum  well  lasers  with 
x=0.3  and  0.4  (compressive  strain),  x=0.5  (lattice  matched),  and  x=0.6  (tensile 
strain)  are  investigated  theoretically  and  experimentally.  In  the  section  3.2, 
the  relative  transition  strength  is  derived  and  the  polarization  anisotropy  is 
discussed  in  the  section  3.3.  The  experimental  results  are  given  in  the  sec- 
tion 3.4.  and  discussion  is  given  in  section  3.5. 

3.2  Transition  Matrix  Element 

For  quantum  well  lasers,  the  carrier  transition  occurs  mostly  near  the 
band  edge  where  the  in-plane  k vector  is  negligible.  Hence  by  estimating  the 
transition  matrix  element  near  the  band  edge,  the  polarization  anisotropy 
can  be  estimated  queilitatively  as  well  as  quantitatively.  The  transition 
matrix  element,  or  the  transition  strength  which  will  be  discussed  later  in 
this  section,  can  be  established  by  defining  the  Bloch  fimctions  in  each  band. 
In  fact,  each  energy  band  can  be  viewed  as  originating  from  the  discrete 
atomic  energy  levels  of  the  isolated  atoms  that  composed  the  crystal  [Kane 
57].  In  this  sense,  the  conduction  band  can  be  thought  as  a remnant  of  an  s 
orbital,  while  the  three  valence  bands  are  the  remnants  of  the  three  p atomic 
orbitals:  p^,  Py,  and  p^-  The  corresponding  Bloch  functions  for  these  orbitals 
are  denoted  as  Ug  (s  orbital),  and  orbital)  where  i 
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denotes  either  x,  y,  or  z.  The  conduction  band  Bloch  function  Ug  has  even 
symmetry  in  all  three  directions  while  Uj-  has  odd  symmetry  along  i direction 
but  even  symmetry  in  the  other  two  directions  so  (u  \u.)  = 0.  However, 

S I 

n 

{u  \p  \u)  is  in  general  non  zero  because  p = (B/i)  (d/dx)  and  hence,  operat- 

S I I 

ing  on  Mj  with  inverts  the  symmetry.  From  these  arguments,  some  sym- 
metry  relations  can  be  written  as  follows  [Zory  93]: 

(upi\u.)  = 0 for  ijtj 

(«^IpK>  = {Ug\p^\u^)  = M (3.1) 

(mJpIm;)  = 0 


where  Mj,  and  ui  indicate  spin-up  and  spin-down  functions.  The  constant  M 
is  defined  here  as  the  basis  fimction  momentum  matrix  element. 

The  valence  band  Bloch  functions  w/j^jfheavy  hole),  «^;j(light  hole),  and 
«so(spin-orbit)  can  be  expressed  in  terms  of  the  linear  combination  of  the 
basis  fimctions  u^,  Uy,  and  [Kane  57,  Lutt  55,  Zory  93].  Defining  the  elec- 
tron’s k vector  to  be  directed  along  z,  the  normalized  valence  band  Bloch  func- 
tions  can  be  written  as  [Zory  93] 


^hh  = 
“/A  ^ 


1 


J2 


^ u +iu  -2u 

X y z 


s 

^/3 


V 


Ux  + lUy  + 


^Ih  = 


^so  = 


u -\u  +2u 

X y z 


(3.2) 


Using  the  above  expressions,  the  matrix  element  which  determines  the 
transition  probability  between  Ug  and  the  basis  functions  (Mj)  can  be  derived. 
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The  conduction  band  effective  mass  along  the  electron’s  k vector  (assume  z 
direction)  is  expressed  as  [Kane  57,  Made  78] 


(3.3) 


where  the  siunmation  sums  over  all  n energy  bands,  and  the  and  are 
the  band  edge  energies  of  each  bands.  By  neglecting  all  but  the  three 
valence  bands  in  the  svunmation,  an  approximate  expression  for  the  conduc- 
tion band  effective  mass  is  obtained  using  the  relation  given  in  Eq.  (3.1)  and 
Eq.  (3.2).  From  Eq.  (3.2),  it  is  noted  that  the  heavy  hole  Bloch  function  does 
not  contain  hence  the  contribution  to  the  sum  is  zero.  Thus,  slim- 

ming over  the  light  hole  and  spin-orbit  band  energy  separations  in  Eq.  (3.3) 

I 

yields 


m* 


J_ 

% 


2IMIV2  _1_  1 

“o  *^3  Eg  3 


1 


(Eg  + A) 


(3.4) 


Since  the  higher  energy  bands  are  neglected  in  this  approximation,  the 
effective  mass  m*  is  expected  to  be  smaller  than  the  true  one  as  illustrated  in 
Fig.3.1.  The  true  conduction  band  effective  mass  can  be  measured  experi- 
mentally to  a good  degree  of  precision  using  cyclotron  resonance  technique. 
By  assiuning  that  m*  is  close  to  the  true  effective  mass  m^,  Eq.  (3.4)  can  be 
rearranged  to  obtain 


^m* 


21  E„  + ^A 

g 3 


(3.5) 
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Figure  3.1  TVpical  III-V  semiconductor  band  structure  schematic. 

The  dashed  line  of  conduction  band  shows  the  approxi- 
mated effective  mass. 
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Although  this  approximation  provides  comprehensive  relationship 
between  the  matrix  element  and  the  energy  gaps,  it  should  noted  that  the 
assumption,  m*  is  close  to  m^.,  may  not  valid  when  the  contribution  from  the 
higher  energy  bands  in  Eq.  (3.3)  is  significant.  For  instance,  in  GaAs,  m*  = 
0.057mQ,  compared  with  the  true  effective  mass  0.067niQ.  However,  it  is  a 
useful  formula  for  a material  that  have  not  been  fully  characterized.  The 
most  accurate  estimate  of  the  transition  matrix  element  can  be  obtained  by 
electron  spin  resonance  techniques. 

Having  obtained  the  magnitude  of  IMI^,  the  transition  matrix  ele- 
ment I M-p  I ^ which  determines  the  transition  probability  between  conduction 
band  and  valence  band  Bloch  functions  can  be  quantified.  The  transition 
matrix  element  I I ^ is  defined  by  [Zory  93] 

|m/  = |<u  Je  ■ p|u,>|^|(F|,|F^>|^  (3.6) 

where  e is  the  unit  polarization  vector,  p is  the  momentum  operator,  and 
Fh(Fg)  is  the  envelope  fimction  of  the  valence  (conduction)  band.  As  shown 
in  Fig.  3.2,  there  are  four  possible  transitions  between  the  spin-degenerated 
conduction  and  valence  bands.  Siunming  over  these  transitions,  Eq.  (3.6) 
becomes 


Uy  UeU„  U, 


(3.7) 


The  envelope  function  overlap  integral  |(F|j|Fi)  is  assumed  equal  to 


unity,  because  for  the  moment  transitions  between  two  bulk  plane  wave 
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Figure  3.2  Diagram  illustrating  the  four  possible  transitions  between 

the  spin-degeneracy  conduction  and  valence  band.  U(.  (uy) 
is  the  spin-up  conduction  (valence)  band  Bloch  function 
and  Uc  (uy)  is  the  spin-down  conduction  (valence)  band 
Bloch  function. 
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electron  states  are  considered.  Later  in  this  section,  more  general  form  of 
the  above  equation  in  the  quantum  well  will  be  discussed.  To  simplify  Eq. 
(3.7),  the  dot  product  between  the  emit  polarization  vector  and  the  electron 
momentum  operator  is  expended  in  the  form  of  Then  by 

using  Eq.  (3.1)  and  Eq.  (3.2),  the  expression  for  is  reduced  to  a very 

simple  form  as  follows: 


M, 


hh 


M, 


so 


= l|M|{|-e^-ie/  + 0 + 0 + 
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Since  +6^  =1,  every  term  of  can  be  replaced  by  1 - 

The  direction  of  k vector  is  again  assmned  z direction.  By  this  assmnption, 

can  be  replaced  by  k • e where  k is  imit  vector  directed  along  k.  Then  Eq. 

2 2 

(3.7)  can  be  rewrited  in  terms  of  M™  /|M|  , or  relative  transition  strength 


2,^2 


as  follows: 
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(3.9) 


Note  that  Eq.  (3.9)  does  not  include  any  reference  frame,  but  the  rela- 
tive orientation  between  the  field  polarization  e and  the  electron  k vector. 
In  Fig.  3.3  and  Fig.  3.4,  the  dependence  of  the  relative  transition  strength 


Transition  Strength 
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Figure  3.3 


Dependence  of  the  transition  strength  on  angle  between 
the  electron’s  k vector  and  the  incident  electric  field  vector 
E for  conduction-heavy  hole  transition. 
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Angle  between  k and  E field 


Figure  3.4 


Dependence  of  the  transition  strength  on  angle  between  the 
electron’s  k vector  and  the  incident  electric  field  vector  E for 
conduction-light  hole  transition. 
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L Z 

|M^|  /lM|  on  angle  between  the  electron’s  k vector  and  the  incident  electric 
field  vector  E is  shown.  For  the  c-hh  transition,  the  relative  transition 
strength  is  maximum  for  the  incident  electric  field  (E)  perpendicular  to  the 
electron  k vector(k)  and  minimmn  for  E parallel  to  k.  The  opposite  holds  for 
c-lh  transition  but  it  should  be  noted  that  the  minimum  transition  strength 


does  not  vanish. 


3.3  Polarization  Anisotropy 

In  this  section,  the  polarization  anisotropy  of  the  quantiun  well  laser  is 
considered  based  on  the  transition  matrix  element  discussed  in  the  previous 
section.  For  QW  structure,  the  envelope  functions  are  constructed  from  two 
or  more  plane  wave  states.  Hence,  cross  terms  between  the  various  plane 
waves  are  contained  which  make  up  the  confined  state  appears  in  Eq.  (3.7) 
for  QW  structure  [Zory  93].  However,  this  complexity  can  he  simplified  by 
restricting  the  region  of  interest  into  near  the  band  edge  and  neglect  the  cross 
terms.  In  this  situation,  the  electron  k vectors  are  quantized  along  certain 
direction.  It  is  illustrated  in  Fig.  3.5  that  how  quantum  confinement  in  a QW 
serve  to  polarized  the  momentum  of  the  band  edge  electrons  along  certain 
directions.  For  QW  lasers  with  heavy  hole  as  a ground  state,  TE  mode  oscil- 
lation prevails  because  the  transition  strength  is  great  when  the  light  is 
polarized  in  the  plane  of  the  QW  (E  parallel  to  the  QW  layer)  as  shown  in  Fig. 
3.5(a).  Conversely,  when  the  light  hole  is  a ground  state,  TM  mode  oscillation 
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Polarization  for  maximum 
transition  strength 


(a) 


Polarization  for  maximum 
transition  strength 


Figure  3.5  Illustration  of  how  quantmn  confinement  serves  to  polarize 

the  momentum  of  band  edge  electrons  for  (a)  C-HH  and  (b) 
C-LH  transitions. 
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is  predominant  because  the  transition  strength  is  maximum  for  E perpendic- 
ular to  the  QW  layer  as  shown  in  Fig.  3.5(b).  The  relative  transition  strength, 


M. 


' I I ^ ^ 

/|M|  , is  calciolated  as  a function  of  the  in-plane  k vector  for  80A 


Ga^Inj^.^P/  (AIq  gGao,4)o.5lno,5P  quantum  well  lasers  with  x=0.3,  0.4,  0.5,  and 


0.6  and  illustrated  in  Fig.  3.6  for  TE  and  Fig.  3.7  for  TM  polarization.  For  TE 
polarization,  the  transition  strength  of  QW  lasers  with  x=0.3,  0.4,  and  0.5  are 
very  strong  near  the  band  edge  because  the  valence  band  ground  state  is 
heavy  hole  state  for  these  devices.  However,  for  TM  polarization,  the  transi- 
tion strength  of  QW  lasers  with  x=0.6  is  much  stronger  than  other  materials 
because  the  light  hole  is  the  valence  band  ground  state. 


3.4  Experiment 

o 

In  this  section,  the  in-plane  spontaneous  emissions  of  80A  GaxIn^.xP/ 
(Alo.6Gao.4)o.5lno.5P  quantum  well  lasers  with  x=0.3  and  0.4  (compressive 
strain),  x=0.5  (lattice  matched),  and  x=0.6  (tensile  strain)  are  investigated. 
Samples  were  grown  by  organometallic  vapor  phase  epitaxy  (OMVPE)  pro- 
vided by  Xerox  Palo  Alto  Research  Center.  The  QW  is  contained  within  a 
typical  diode  laser  structure  as  shown  in  Fig.  2.4.  Confining  layers  are 
undoped  (AIq  gGao  ^lo.sIno.sP,  while  p-  and  n-cladding  layers  are  AIq  glng  5P. 
Diodes  were  100pm  wide  by  500pm  long  stripes  with  cleaved  facets,  which 
were  probe  contacted  and  pulse  tested. 

The  emitted  light  was  collected,  and  focused  onto  the  entrance  slit  of  a 
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Figure  3,6  Relative  transition  strength  of  conduction  to  first  valence 

band  TE  polarization  for  80A  Ga^Inj^.^P/ 

(AIq  6Grao.4)o.5^^0.5P  quantiun  well  lasers  with  x=0,3,  0.4, 
0.5,  and  0.6 


Relative  Transition  Strength 


50 


Figure  3.7  Relative  transition  strength  of  conduction  to  first  valence 

band  TM  polarization  for  80A  GaxIn^.xP/ 
(Alo.6Gao.4)o.5^^^0.5P  quantum  well  lasers  with  x=0.3,  0.4, 
0.5,  and  0.6 
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spectrometer,  using  a microscope  objective.  A Glan-Thompson  Air-Spacing 
Polarizer  (clear  aperture  10mm)  was  used  to  analyze  the  spontaneous  emis- 
sion in  the  desired  polarization.  Shown  in  Fig.  3.8,  Fig.  3.9,  and  Fig.  3.10  are 
the  TE  mode  and  TM  mode  electroluminescence  (EL)  spectra  of  the  three 
samples  for  x=0.4,  0.5,  and  0.6,  taken  at  a current  density  of  about  200,  130, 
and  100mA  respectively.  As  the  GaP  alloy  content  in  the  QW  is  increased 
from  x=0.4  to  x=0.6,  the  band  gap  increases  and  there  is  a net  shift  to  shorter 
wavelength.  The  polarization  of  the  EL  spectra  is  anisotropic,  with  the  e-hh 
transition  TE  polarized,  and  e-lh  transition  predominantly  TM  polarized. 

For  the  compressively  strained  QWs  (x=0.3  and  x=0.4),  the  TE  mode 
emission  is  strongest,  because  the  hh  band  is  ground  state.  The  opposite  is 
true  for  tensile  strain  (x=0.6),  in  which  case  the  TM  mode  emission  is  most 
intense  because  the  valence  band  edges  have  become  reversed.  At  very  low 
current  densities  Gess  than  5A/cm^,  for  which  b£md  filling  is  negligible  and 
the  Boltzmann  approximation  to  the  Fermi  distribution  applies),  the  relative 
intensity  between  the  TE  and  TM  mode  emission  (I-j-e/Itm)  follows  the  ratio 
between  the  occupancy  factors  for  the  Ih  and  hh  bands,  i.e.,  exp(AE'rM-TE/^'I’)> 
where  AE^m-te  is  the  hh-lh  band  edge  difference,  and  kT  is  the  thermal 
energy;  this  is  summarized  in  Table  3.1.  The  case  near  lattice  match  (x=0.5) 
has  the  Ih  and  hh  band  edges  nearly  degenerate,  resulting  in  similar 
wavelength  and  intensity.  Exact  degeneracy  of  the  Ih  and  hh  band  edges 
requires  slight  biaxial  tension,  however.  For  this  QW  thickness,  the 
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Wavelength  (nm)  • 


Figure  3.8  Electroluminescence  of  80A  CaxIn^.^P/ 

(Alo.eGrao.4)o.5lno.5P  quantum  well  lasers  (x=0.4)  analyzed 
in  both  TE  and  TM  polarizations,  at  a current  I=200mA. 
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Wavelength  (nm)- 


Figure  3.9  Electroluminescence  of  80A  GaxIn^.xP/ 

(AIq  6Gao.4)o.5lno.5P  quantum  well  lasers  (x=  0.5)  analyzed 
in  both  TE  and  TM  polarizations,  at  a current  1=  130mA. 
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o 

Figure  3.10  Electroluminescence  of  80A  GuxIni.xP/ 

(Alo.eGao .4)0.51110.5?  quantum  well  lasers  (x=  0.6)  ana- 
lyzed in  both  TE  and  TM  polarizations,  at  a current 
I=100mA. 
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Table  3.1 

The  comparision  of  the  TE  and  TM  emission 


X 

^TE/TM 

AEtm-TE 

AEtm-te 

(Experiment) 

(Theory) 

0.3 

6.5 

42.0  meV 

32.1meV 

5.1 

0.4 

4.2 

33.2  meV 

32.7meV 

3.6 

0.5 

0.95 

-1.0  me  V 

6.1  meV 

» 

0.96 

0.6 

0.56 

-12.4  meV 

-30.1meV 

m 

0.62 

56 


appropriate  composition  is  x=0.53  [Bour  93].  The  energy  separation  between 

the  TE  and  TM  mode  spontaneous  emission  peaks  reflects  the  energy 

difference  between  the  light  and  heavy  hole  band  edges.  Contributing  to  this 

difference  are  both  the  strain  induced  splitting  of  the  Ih  and  hh  bands,  and 

the  different  quantxun  shifts  of  the  two  valence  bands  [Pan  88,  Kond  91, 

Wang  90,  Van  89].  The  measured  energy  separation  (AEtm-te)  well  as  the 

calculated  value  is  given  in  the  Table  3.1.  The  theoretical  value  was 

calculated  based  on  the  c-hh  and  c-lh  n=l  energy  level  difference  calculated  in 

Chapter  II.  Since  the  origin  of  the  TE-TM  peak  wavelength  separation  is  the 

hh  and  Ih  splitting  at  the  band  edge  due  to  the  shear  strain  deformation,  the 

« 

polarization  anisotropy  can  be  predicted  simply  by  calculating  the  energy 
split  defined  in  Eq.  (2.12).  There  was  significant  disagreement  between  the 
theory  and  the  experiment  for  x=0.3  and  0.4  and  this  will  be  discussed  in 
detail  in  the  next  section. 

The  current  density  dependence  of  the  TE  and  TM  mode  spontaneous 
emission  peak  wavelength  are  shown  in  Fig.  3.11.  The  wavelength 
separation  between  TE  and  TM  mode  emission  is  constant  at  low  current 
densities.  At  elevated  current  densities,  however,  bandfilling  begins  to  occur, 
and  the  wavelength  of  the  fundamental  transition  shortens  with  increasing 
injection  current.  This  bandfllling  effect  is  most  pronoimced  in  the  devices 
with  x=0.3  and  x=0.4,  suggesting  that  the  in-plane  effective  mass  of  heavy 
holes  is  reduced  in  such  QWs  under  biaxial  compression.  Furthermore,  the 
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Figure  3. 11  Current  density  dependence  of  spontaneous  emission  peak 

wavelength  for  80A  CaxIn^.^P/  (Alo  eGaQ  ^lo.sIno  5P  quan- 
tum well  lasers  (x=0.3, 0.4,  0.5,  0.6)  analyzed  in  both  TE 
and  TM  polarizations. 
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threshold  current  density  is  also  lowest  in  these  devices,  because  population 
inversion  occurs  at  a lower  QW  carrier  density  when  the  in-plane  effective 
mass  of  the  hh  band  is  reduced. 

For  the  devices  with  x=0.5  and  0.6,  a shift  to  shorter  wavelengths  at 
high  current  density  is  not  apparent.  This  suggests  that  the  in-plane  hole 
effective  mass  might  be  greater  than  for  the  compressively  strained  QWs. 
Thresholds  are  also  higher  in  these  devices,  increasing  to  400A/cm^  for  the 
TM  mode  device  with  x=0.6.  Presumably,  this  is  due  to  both  their  higher  in- 
plane effective  masses,  and  the  diminishing  electron  confinement  at  short 
wavelengths.  This  will  be  more  clarified  in  the  next  section.  Also  in  Fig. 
3.11,  the  reversal  in  the  hh  and  Ih  band  edges  is  evident,  as  the  biaxial  QW 
strain  progresses  from  compression  to  tension. 

Above  threshold,  the  polarization  of  stimulated  emission  was  predomi- 
nantly the  same  as  that  of  the  ftmdamental  transition.  The  spectra  for 
x=0.4,  and  0.6  devices  at  the  threshold  and  high  injection  for  the  alternative 
polarization  (TM  for  x=0.4  and  TE  for  x=0.6)  are  shown  inFig.  3.12,  and  Fig. 
3.13.  The  lasing  peak  wavelength  appeared  at  the  same  wavelength  with 
the  ftmdamental  transition  lasing  peak  rather  than  around  the  middle  of  the 
spontaneous  emission  spectrum  and  hence  the  polarization  state  of  laser 
oscillation  was  not  pure. 


Intensity  (a.u) 


59 


Wavelength  (nm) 


Figure  3.12  The  measured  TM  spectrum  of  Gao  ^Ino  gP  80A  QW  lasers 

near  the  threshold  (300mA)  and  well  above  the  threshold 
(350mA) 


Intensity  (a.u) 
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Wavelength  (nm) 


Figure  3.13  The  measured  TE  spectrum  of  Gao.6l1to.4P  QW  lasers 

near  the  threshold  (800mA)  and  well  above  the  threshold 
(900mA) 
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3.5  Discussion 

It  has  become  apparent  that  the  lowering  of  valence  band  effective 
mass  can  have  considerable  benefits  for  semiconductor  lasers  [Yabl  86].  The 
intentional  incorporation  of  strain  in  QW  layer  has  become  popular  technique 
for  reduction  of  the  valence  band  effective  mass.  The  compressive  strain  as 
well  as  tensile  strain  causes  a considerable  reduction  of  the  in-pleme  effective 
mass  [Krij  92].  As  a result,  the  hole  density  states  is  reduced  and,  in  the 
ideal  case,  becomes  comparable  to  the  electron  density  of  states.  In  such  con- 
dition,  the  carrier  density  required  for  degeneracy  would  be  minimized  and 
the  Bernard  -Duraffourg  condition  [Bern  61]  (Fq  > hv  > Eg  cotJd  be  sat- 
isfied at  a low  carrier  density  with  accompan}dng  benefits  in  terms  of  thresh- 
old current  density  [Yabl  88]. 

The  energy  vs.  in-plane  k dispersion  under  the  strain  effect  can  be  cal- 
culated by  estimating  the  deformation  potentials  in  the  QW  layer  and  insert- 
ing into  the  Hamiltonian  as  discussed  in  chapter  II.  The  results  were  shown 
in  Fig.  2.10  and  Fig.  2.11.  By  assuming  that  the  peak  emission  of  light  came 
from  the  band  edge  transition  from  conduction  electron  to  the  heavy  hole  (TE) 
and  light  hole  (TM),  the  peak  TE  and  TM  wavelength  difference  (AEtm-te) 
were  predicted  as  listed  in  Table  3.2.  As  discussed  in  the  previous  section, 
however,  there  was  considerable  disagreement  in  the  AE'pi^.'j'E  between  the 
experiment  and  the  theory.  In  determining  AE-pM-TE>  most  important 
factor  is  the  shear  deformation  potential  because  the  hh  and  the  Ih  state  are 
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split  by  this  strain.  The  splitting  (S)  is  defined  as  [Poll  68,  Gavi  70,  Chan  77] 


S = AEjj  (shear)  -AEj  (shear) 


= 6E 


1 (5Eooi) 


(3.10) 


001  2 
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where  SEqqj  is  the  linear  splitting  of  the  naultiplet  for  stress  parallel  to 


[001]  direction  defined  as  6E 


001 


= 2b(e  -e  ) and  b is  the  shear  deforma- 


tion  potential.  The  strain  tensors  e and  e are  given  as 
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Combining  Eq.  (3.10)  and  Eq.  (3.11); 
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(3.12) 


The  second  term  in  Eq.  (3.12)  is  often  neglected  for  large  A^  such  as  GaAs 
(Ao«340meV),  but  must  be  included  for  small  A^  such  as  GalnP  (Ao~110meV). 

a 

The  parameters  used  in  Eq.  (3.12)  for  GaxIn^.xP  were  linearly  interpolated 

from  host  alloys  GaP  and  InP.  The  only  factor  in  Eq.  (3.12)  which  depends 
strongly  on  x is  Aa/aq^  where  Aa  is  the  lattice  constant  difference  between  the 


quantum  well  and  the  barrier  layer  Ubarrigr  ‘ ^qw*  Eig.  3.14.,  abarrier/^ 


^qw 


VS.  X is  plotted.  The  solid  line  is  derived  from  linear  interpolation  of  the  lat- 
tice constants  of  GaP  and  InP  and  the  triangles  are  determined  so  that  the 
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Ga  mole  fraction  x 


Figure  3.14  The  ratio  of  abamer^^w  of  GaxIni.xP  as  a function  of  Gallium 

o 

mole  fraction  x assuming  abarrier=5.6533A.  The  solid  line  is 
linearly  interpolated  data  from  GaP  and  InP  and  the  triangles 
are  determined  so  thst  the  calculated  splitting  S 3uelds  the 
same  AE^e-tm  experiment. 
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calcialated  splitting  3delds  the  same  AE-te-tm  experiment.  For  x=0.53, 

QW  layer  is  nearly  lattice  matched  to  the  barrier  layer.  As  can  be  seen  in 

« 

Fig.  3.14,  the  ratios  (abarrie/^w)  modified  aq^^  for  x=0.4  (about  0.8%  lat- 
tice mismatch)  and  x=0.6  (about  0.6%  lattice  mismatch)  are  close  to  the  inter- 
polated data.  However,  for  x=0.3  (about  1.6%  lattice  mismatch),  there  is  a 
significant  deviation  from  the  interpolated  data.  This  implies  that  the  strain 
acting  on  the  QW  layer  may  not  be  simply  derived  from  Eq.  (3.12)  but  from  a 
more  complicated  form  for  lattice  mismatches  beyond  certain  range.  As 
shown  in  chapter  II,  the  energy  shift  in  a strained  QW  was  defined  based  on 
the  recoil  force  due  to  the  lattice  mismatch  by  analogy  with  Hook’s  law. 
Although  a strained  QW  layer  can  avoid  the  development  of  interfacial  dislo- 
cations if  the  thickness  is  less  than  a certain  critical  value,  the  strain  effect 

« 

pan  be  more  complicated  for  a highly  strained  layer.  For  the  1.6%  lattice 
mismatch  case  (x=0.3),  it  may  be  that  Eq.  (3.12)  is  no  longer  valid.  On  the 
other  hand,  it  may  that  the  linear  interpolation  shown  in  Fig.  3.14  for  GalnP 
is  incorrect.  In  any  event,  polarization  anisotropy  spectroscopy  provides  a 
means  for  studying  these  new  strained  QW  structure. 


CHAPTER  IV 

TEMPERATURE  DEPENDENCE  OF  THE  THRESHOLD  CURRENT 
DENSITY  OF  CdZnSe  BLUE-GREEN  QUANTUM  WELL  LASERS 

4.1  Introduction 

Blue  semiconductor  lasers  have  great  advantages  over  commonly  used 
infrared  (IR)  lasers.  For  example,  optical  readout  and  storage  densities  can 
be  increased  by  a factor  of  three  if  one  substitutes  a blue  laser  for  an  IR  laser. 
Despite  the  advantages,  the  development  of  ZnSe-based  blue  semiconductor 
lasers  has  been  slow  because  of  the  difficulty  in  growing  good  p-type  material. 
Recently,  a new  p-type  growth  technique  using  a Nitrogen  plasma  during 
molecular  beam  epitaxy  (MBE)  has  been  developed  [Park  90].  This 
technique  was  used  to  make  the  first  blue-green  semiconductor  injection  laser 
by  researchers  from  3M  [Haas  91].  Since  then,  various  kinds  of  II- VI 
material  have  been  investigated  for  better  performance  and  reliability. 
These  achievements  give  us  confidence  that  we  are  on  the  right  path  toward 
compact,  room  temperature  operating  blue  diode  lasers.  Recently,  Okuyama 
et  al  [Okyu  91]  proposed  ZnMgSSe  as  a new  cladding  layer  material.  Gaines 
et  al  [Gain  93]  proved  its  superiority  by  achieving  the  lowest  room 
temperature  threshold  current  density  and  the  highest  operating 
temperature  using  a Cdo.2Zno.8Se  quantum  well,  ZnSo.o6Seo.94  barriers,  and 
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Zhq  gMgo  iSo.iSeo  g-  cladding  layers.  Using  this  same  QW  structure,  Haase 
et  al  [Haas  93]  from  3M  reported  more  recently  the  record  low  threshold 
current  with  a buried  ridge  waveguide  design.  In  this  chapter,  we  calculate 
the  threshold  current  density  of  CdZnSe/ZnSe  strained  single  quantmn  well 
diode  lasers  over  a wide  temperature  range  and  compare  with  the  published 
experimental  data. 


4.2  Spontaneous  emission  rate  and  gain  spectra 


Using  the  E vs.  k dispersion  derived  in  the  chapter  2,  the  quasi-Fermi 
energy  levels  were  calculated.  The  carrier  densities  of  the  conduction  and 
valence  band  can  be  written  as 


m kpT 

N = -^Sln{l  + exp[-(E^„-Ef^)/(kBT)]} 
TcR  L n 

P = 2EJ!""p“(V  n - f V [e.„  (k,y)  ] } 

n ^ 


(4.1) 


where  Efg  (Efy)  is  the  quasi-Fermi  level  of  the  conduction  (valence)  band,  L is 
the  thickness  of  the  QW,  fy  is  the  fermi  distribution  function  for  the  valence 
band,  fy=l/[l-i-exp{E  - Efyl/kgT],  and  Egjj  (Eyj^)  are  the  quantized  energy  level 
of  the  conduction  (valence)  band.  The  density  of  the  state,  p (k  ) , for  two 

xy 

kxv  1 

dimensional  electron  system  is  given  by  p (k  ) = • - . The  summation 

in  Eq.  (4.1)  should  be  done  over  all  quantized  energy  states,  even  on  the 
unboimded  states  which  have  higher  energies  than  QW  barrier  as  pointed  out 
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by  Nagarajan  et  al  [Naga  89],  for  a numerical  accuracy.  This  is  important 
especially  for  high  injection  case.  The  quasi  Fermi  levels  with  respect  to  the 
carrier  density,  N or  P,  were  found  by  assuming  the  charge  neutrality,  N=P, 
and  the  results  are  shown  in  Fig.  4.1  and  Fig.  4.2. 


The  spontaneous  emission  rate  (rgp)  is  given  by  [Zory  93] 


^sp(hv)  = 


2tc 


^eV 


If 


R Amo>' 


2R  (co/e)  ^ 
4o)  J 


M-rl  (4.2) 


where  e is  the  electron  charge,  m^  is  the  free  electron  rest  mass,  £ is  the  elec- 

2 

trie  permittivity,  (hv)  is  the  optical  mode  density,  |M.j,|  is  the  transition 
matrix  element  defined  in  Eq.  (3.7),  f^Cfy)  is  the  fermi  distribution  of  conduc- 
tion (valence)  band,  and  p . (hv)  is  the  reduced  density  of  states  given  by 


Pfed  (R^) 


.dEeh/dk^ 


(4.3) 


E = hv 


Since  one  cannot  use  the  closed  form  expression  for  p ^^(hv)  for  the 
nonparabolic  valence  band  case,  we  calculated  numerically  from  the  energy 
vs.  in-plane  k vector  dispersion  curve. 

The  spectral  gain  is  given  by  [Zory  93] 


g(hv) 


27ce^R 

eQ^ncm^hv 


MlfpredC"') 


(4.4) 


where  n is  a refractive  index  parameter,  c is  the  speed  of  light,  and  X is  the 
wavelength.  The  calculated  spontaneous  emission  rate  and  the  gain  spectra 
for  Cdo.2Zno,8Se  are  shown  inFig.  4.3,  and  Fig.  4.4. 
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Figure  4.1  The  quasi-Fermi  energy  levels  for  Cdo.2Zno.8Se  65A  QW 

structure.  The  zero  energy  level  is  the  conduction  band 
edge  of  the  QW. 
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Figure  4.2  The  quasi-Fermi  energy  levels  for  Cdo  2Zno  gSe  65A  QW 

structure.  The  zero  energy  level  is  the  valence  band  edge 
of  the  QW. 
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Figure  4.3  Spontaneous  spectra  of  CdZnSe  65A  QW  diode  lasers. 
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Figure  4.4 


Spectral  gain  of  CdZnSe  65A  QW  diode  laser 
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4.3  Intraband  Relaxation  and  Line  Shape  Function 

Line  shape  functions,  representing  homogeneous  broadening  due  to 
carrier  scattering,  are  defined  as  response  functions  of  dipole  phase  damping 
[Yama  87].  Line  shapes  for  optical  transitions  play  an  important  role  in  esti- 
mating optical  properties  such  as  a gain  and  a spontaneous  emission  spectra. 
It  is  widely  accepted  that  the  line  shape  function  for  semiconductors  is  the 
Lorentzian  fimction  [Yama  81],  in  analogy  with  the  gas  laser  theories.  How- 
ever, when  the  Lorentzian  fimction  is  applied  to  QW  semiconductor  lasers, 
the  gain  coefficients  are  underestimated  and  an  unnatural  absorption  region 
appears  at  photon  energies  below  the  band  gap  in  gain  spectra  [Lim  91]. 
Recently,  line  shape  functions  for  QW  structures  were  theoretically  analyzed 
by  taking  non-Markovian  relaxation  processes  into  accoimt  [Yama  87,  Asad 
89,  Ohto  91].  It  has  been  foimd  that  the  line  shape  fimctions  in  low  dimen- 
sional  systems  have  a strong  convergent  characteristic  and  asymmetry. 
Since  the  scattering  time  of  ZnSe-based  material  is  known  to  be  very  short 
compared  to  conventional  III-V  materials  such  as  GaAs  [Nurm  93],  the 
proper  choice  of  the  line  shape  function  is  very  important. 

The  line  shape  fimction  L(hv)  is  given  by  the  Fourier  transform  of  the 
optical  dipole  moment  P(t)  [Yama  87]: 

L (hv)  = ^ (t)  (4.5) 


The  dipole  moment  P(t),  which  is  damped  due  to  carrier  scattering,  can  be 
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written  as 


TiCOnt  -Fft) 

P(t)  = P(0)e  e 


(4.6) 


where  cOq  is  the  resonant  angular  frequency  and  F(t)  is  the  phase  damping 


factor.  Here,  it  is  important  to  include  the  non-Markovian  processes  which 
bring  about  the  nonlinear  time  dependence  of  the  damping  factor  F(t).  This 
gives  rise  to  the  crucial  difference  between  the  present  line  shape  and  the 
conventional  Lorentzian  one.  Without  tracing  the  temporal  behavior  of 
dipole  momentum  P(t),  Asada  obtained  the  Fourier  transform  of  P(t)  using 
the  temperature  Green’s  method  [Asad  89].  In  that  work,  hole-hole,  hole- 
electron,  and  carrier-longitudinal  optical  (LO)  phonon  scattering  time  were 
analyzed.  It  has  been  foimd  that  hole-hole  scattering  is  dominant  among  the 
various  scattering  due  to  its  heavy  effective  mass.  Moreover,  an  approximate 
form  of  the  line  shape  L(hv)  for  numerical  calculation  has  been  derived  for 
hole-hole  scattering  (which  dominates  the  various  scattering  processes)  as  fol- 
lows [Zory  93]: 
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where  x^  is  the  hole-hole  scattering  time,  Efy  is  the  quasi-Fermi  energy  level, 
Lyj  is  the  effective  quantiim  well  width  for  jth  valence  band  holes,  n2D  is  the 
two  dimensional  carrier  density,  e is  the  electrical  permittivity,  mg(my)  is  the 
electron  (hole)  effective  mass,  Eg  (E^)  is  the  energy  level  of  the  conduction 
(valence)  bemd,  and  ky  is  the  valence  band  in-plane  k vector.  As  defined  in 
Eq.  (4.8),  r (hv  -Epy)  is  a function  of  the  photon  energy.  However,  for  the 
conventional  Lorentzian  line  shape  function,  F (hv  - E ) is  simply  R/Xj^^ : 


L(hv-E^„)  = 1 


R/X: 


in 


it(hv-E^,)%rE/Xj„ 


(4.9) 


where  x-  is  the  intraband  relaxation  time.  The  comparison  between  the 
Lorentzian  and  modified  Lorentzian  line  shape  function  for  CdZnSe  65A  QW 


is  illustrated  in  Fig.  4.5.  As  can  be  seen  in  the  figure,  the  modified 


Line  Shape  L(hv-Ecv)  (eV) 
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Figure  4.5  Calculated  results  of  Lorentzian  and  modified  Lorentzian 

o 

line  shape  function  at  300K  for  65A  CdZnSe  quantum  well 
lasers.  The  scattering  time  is  assumed  60  femtoseconds. 
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Lorentzian  line  shape  is  narrower  and  the  low  energy  tail  is  much  steeper 

than  the  Lorentzian;  it  decays  exponentially.  In  calculating  the  line  shapes 

in  Fig.  4.5,  the  scattering  time  of  60  femtoseconds  is  assumed  from  a recent 

experimental  result  [Nurm  93].  This  is  a significantly  shorter  scattering 

time  compared  to  conventional  III-V  materials.  Since  the  scattering  time  is 

a function  of  temperature,  the  line  shape  changes  at  different  temperatures. 

The  variation  of  the  scattering  time  with  temperature  is  caused  by  the 

thermal  distribution  of  the  carriers  and  the  inverse  screening  length  a which 

is  defined  in  Eq.  (4.8).  For  the  calculation  of  the  scattering  time  of 

Cdo  ^Zno  gSe  QW  lasers  as  a function  of  temperature,  the  data  point  of  60 

femtoseconds  at  300K  was  assumed  and  extrapolated.  The  results  are 

shown  in  Fig.  4.6.  The  lineshapes  for  various  temperatures  were  then 

calculated  for  hole-hole  scattering.  The  parameters  used  for  Cdo  ^Zno  ^Se 

are  listed  in  Table  4.1  and  the  lineshape  for  n=l  transition  at  85K,  150K,  and 

300K  are  shown  in  Fig.  4.7.  The  line  shape  fimctions  are  normalized  so  that 

« 

the  area  under  the  curve  is  unity.  As  could  be  expected,  the  line  shape 
fimction  becomes  narrower  at  low  temperature  which  implies  less  spectral 
broadening  of  optical  gain  and  spontaneous  emission  rate. 

4.4  Temperature  Dependence  of  The  Threshold  Current  Density  of 

CdZnSe  Blue-Green  Quantum  Well  Lasers 

Using  the  modified  line  shape  function  in  section  4.3,  the  spectral  gain 
was  convolved  in  order  to  take  into  acco\mt  the  intraband  relaxation  which 
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Figure  4.6  The  calculated  temperature  dependence  of  scattering  time  of 

CdZnSe  6SA  QW  diode  lasers.  The  data  point  of  60  femto- 
seconds at  300K  was  used  as  a fitting  parameter. 
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Figure  4.7  The  modufied  line  shape  functions  at  various  tempera- 
tures for  CdZnSe  65A  QW  laser. 
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gives  rise  to  the  spectral  broadening  of  each  transition.  The  intraband  relax- 
ation effect  is  taken  into  account  in  the  gain  calculation  as  follows: 

g,(hv)  =^^g(E„)L(hv-E^,)dE^„  (4.10) 

where  gc(hv)  is  the  spectral  gain  including  the  effect  of  the  intraband 
relaxation,  Eg  is  the  energy  gap  between  the  quantized  energy  levels  of  the 
conduction  and  valence  band,  hv  is  the  photon  energy,  E(.y  is  the  transition 
energy,  and  gfE^y)  is  the  spectral  gain  without  the  intraband  relaxation 
which  is  expressed  in  Eq.  (4.4).  The  convolved  gain  g^.  at  a specific  photon 
energy  hvg  is  determined  by  the  stun  of  the  imconvolved  gain  g mvJtiplied  by 
the  lineshape  function  L.  As  a weighting  factor,  the  line  shape  function  L is 
maximum  at  the  photon  energy  hvo  and  decreasing  for  photon  energies  away 

from  hV().  However,  g(hv)  has  a negative  value  at  large  photon  energies  as 

$* 

can  be  seen  in  Fig.  4.4  and  the  negative  values  are  also  included  in  the 
summation.  Those  negative  terms  are  diverging  because  the  weighting 
factor  L approaches  zero  faster  and  hence,  the  total  sum  of  the  negative  value 
is  finite.  The  contribution  of  the  finite  negative  sum  gives  rise  to  the 
decrease  of  the  peak  gain.  Also  it  should  be  noted  that  the  gain  does  not 
necessarily  vanish  for  the  photon  energies  less  than  E„  which  was  forbidden 
for  unconvolved  gain  g.  The  convolved  and  unconvolved  gain  spectra  of 
Cdo.2Zno.8Se  are  shown  in  Fig.  4.8.  It  would  be  worth-while  to  consider  the 
physical  significance  of  the  numerical  result  of  the  convolution.  Without 
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Figure  4.8  The  convolved  and  xmconvolved  spectral  gain  of  CdZnSe  65A 

QW  diode  lasers  at  300K  for  carrier  density  lel9cm'^. 
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convolution,  the  gain  is  determined  solely  by  the  carrier  distribution  in  the 
band  and  the  transition  matrix  element  which  was  discussed  in  chapter  III. 
Without  scattering,  the  carrier  distribution  in  k space  does  not  change  in  the 
steady  state.  In  this  situation,  the  carrier  transition  for  the  specific  photon 
energy  hvg  only  occurs  between  electrons  and  holes  for  corresponding  in- 
plane  k vector  ko.  However,  with  the  carrier  scattering  included,  the 
electrons  and/or  the  holes  which  were  not  originally  at  kQ,  contribute  to  the 
transition  for  the  photon  energy  hvQ  as  illustrated  in  Fig.  4.9.  Suppose  that 
a transition  with  larger  photon  energy  hv’  at  k’  gives  rise  to  a negative  gain 
without  convolution  as  shown  in  Fig.  4.10.  The  electrons  or  holes  which  were 
not  involved  in  the  transition  at  k=kp  when  the  scattering  was  neglected  now 
contribute  to  the  gain  calculation  through  the  convolution  integral  Eq.  (4.10). 
Some  of  them  may  have  larger  k value  than  k’  -and  hence  correspond  to  a 
larger  photon  energy  than  hv’  which  gives  rise  to  a negative  gain.  These 
negative  gain  terms  are  also  weighted  and  included  in  the  gain  for  the  photon 
energy  hv^.  This  is  why  the  convolved  gain  is  smaller  than  the  unconvolved 
one  in  spite  of  the  momentum  and  the  energy  conservation  of  the  carrier 
scattering  process.  Another  important  feature  in  Fig.  4.8  is  that  near  the 
photon  energy  where  the  gain  is  zero  (E’=hv’  in  Fig.  4.10),  the  convolved  and 
unconvolved  gain  spectra  coincide  as  opposed  to  Chinn  et  al  [Chin  89].  This 
is  a logical  consequence  because  armmd  k’,  the  contribution  from  either  side 
(positive  and  negative  gain)  due  to  the  scattering  is  about  the  same.  If  the 
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Figure  4.9  Schematics  of  carrier  scattering  in  conduction  and 

valence  band. 
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Figure  4.10  Schematic  illustrtions  of  the  effect  of  intraband  scattering  on 

spectral  gain. 


85 


line  shape  function  is  perfectly  symmetric,  the  cross  point  k’  will  be  exactly 
the  same  provided  that  the  gain  decreases  linearly  around  k’.  Although  the 
modified  line  shape  fimction  used  here  is  as3mimetric,  it  converges  much 
faster  than  a conventional  Lorentzian  so  the  cross  points  are  still  quite  close. 
If  the  line  shape  function  did  not  converge  fast,  the  large  negative  part  of 
gain  for  k larger  than  k’  would  exceed  the  positive  contribution.  In  this  case, 
the  cross  point  wordd  have  a smaller  k than  the  unconvolved  one. 

In  order  to  predict  the  threshold  current  density,  the  convolved  gain  gj, 
and  the  corresponding  curreny  density  (J)  are  determined  the  g-J  method  as 
following  procedure: 

1)  Specify  carrier  density  N and  temperature  T. 

2)  Calculate  the  spectral  dependence  of  the  spontaneous  recombination 
rate  (rgp(hv))  and  the  convolved  spectral  gain  (ggfhv))  and  determine 
the  peak  gain  gpeak- 

3)  Integrate  rgp(hv)  to  find  the  current  density  J corresponding  to  that 
value  of  gpeak- 


The  spontaneous  recombination  rate  is  related  to  current  density  J by 


*sp 


(4.11) 


where  d is  the  active  layer  thickness  and  T)  is  the  efficiency  for  converting  J 

sp 

into  spontaneous  emission  in  the  active  layers.  We  assumed  =1  which  is 
a good  approximation  for  long  cavity  laser  at  low  temperature.  The 
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schematic  illustrations  of  the  g-J  method  are  shown  in  Fig.  4.11,  From  these 
relationships,  we  derived  the  peak  gain  vs.  J for  CdZnSe  as  shown  in  Fig. 
4.12 

Based  on  the  g-J  diagram,  we  found  by  estimating  the  threshold 
mode  gain  (G^h)  which  is  given  by 

G.h  = rg„  = a..iln^^  (4.12) 

where  a-  is  the  mode  loss  coefficient,  L is  the  cavity  length,  T is  the  mode 
confinement  factor  in  the  QW,  and  R^,  R2  is  the  facet  refractivity.  The 
calculated  mode  intensity  profile  for  the  CdZnSe  structure  in  Fig.  2.5  is 
shown  in  Fig.  4.13.  Those  parameters  used  in  the  calculation  are  a.  =8.1cm‘^, 
L=1000  pm,  Ri=R2=0.23,  and  F =0.017.  In  our  estimation,  was  18.6cm'^ 
and  hence  g^jj  was  1100  cm'^.  In  terms  of  the  optical  confinement  point  of 
view,  r of  the  device  structure  shown  in  Fig.  2.5  was  enhanced  by  about  25% 
compared  to  the  first  demonstrated  blue-green  laser  (F  =0.0136)  [Haas  91] 
owing  to  the  lower  refractive  index  of  the  ZnSSe  cladding  layer. 

The  threshold  current  density  was  determined  by  the  g-J  curve  and  the 

o 

estimated  g^h-  The  theoretical  values  of  J^h  at  various  temperatures  are 
compared  with  experiment  as  shown  in  Fig,  4.14.  The  predicted  values  agree 
well  with  the  experiment  up  to  300K.  Above  300K,  the  measured  J^^  was 
much  higher  than  predicted,  the  discrepancy  being  due  most  likely  to  carrier 
leakage.  One  striking  result  is  that  the  difference  of  the  contour  between  the 
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Figure  4.11  Schematics  of  the  g-J  method;  (a)  the  spectral  gain  at  vari- 
ous carrier  densities  (b)  the  relation  between  J and  rgp(hv) 
(c)  the  g-J  curve 
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Figure  4.12 


The  peak  gain  vs.  current  density  of  Cdo  ^ZnQ  gSe  QW 
lasers  at  various  temperatures. 
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Figure  4. 13  The  mode  intensity  and  the  refractive  index  profile  of 

Cdo.2Z  .Q.gSe/ZnSSe  6SA  QW  lasers. 
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Figure  4.14  The  calculated  and  measured  temperature  dependence  of 

threshold  current  density  of  Cdo.2Zno.8Se  65A  QW  lasers. 
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convolved  £ind  the  unconvolved  curves;  the  convolved  curve  is  convex  which  is 
consistent  with  experimental  data  while  the  unconvolved  one  is  concave 
[Park  93].  This  is  mainly  due  to  the  temperature  dependence  of  the  carrier 
scattering  time.  As  shown  in  Fig.  4.6,  the  scattering  time  is  exponentially 
decreasing  as  temperature  increases  and  results  in  a wider  line  shape  at 
higher  temperature  as  illustrated  in  Fig.  4.7.  As  a consequence,  the  avail- 
able peak  gain  diminishes  by  more  at  higher  temperatures  and  hence,  more 
current  is  required  in  order  to  reach  the  threshold  gain  which  in  turns 
increase  the  threshold  current  density. 

9 

4.5  Stunmary 

The  conventional  electron  hole  plasma  (EHP)  model,  that  optical  gain 
is  derived  from  k-conserved  stimulated  processes  between  free  electrons  and 
holes  in  quantum  well  subbands,  has  been  used  to  predict  the  threshold  cur- 
rent density  of  CdZnSe  single  quantum  well  lasers  and  good  agreement  with 
experiment  was  obtained  up  to  300K  provided  a hole-hole  scattering  time  of 
60  femtoseconds  is  assumed.  Above  300K,  the  measured  threshold  current 
density  was  much  higher  than  predicted,  the  discrepancy  being  due  most 
likely  to  carrier  leakage.  A modified  Lorentzian  line  shape  function  was  used 
to  include  the  intraband  relaxation  process  caused  by  carrier  scattering.  The 
modified  Lorentzian  line  shape  fimction  for  the  two  dimensional  system  (QW) 
was  much  narrower  than  the  conventional  Lorentzian.  The  line  shape  was  a 
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function  of  temperature  and  showed  a stronger  convergence  characteristic  at 
low  temperatures. 


CHAPTER  V 

DUAL  PEAK  WAVELENGTH  OPERATION  OF  AlGaAs  AND  CdZnSe  QW 

LASERS 

5.1  Introduction 

The  temperature  dependence  of  threshold  current  in  quantmn  well 
lasers  can  be  quite  complex.  For  example,  in  certEiin  temperature  ranges, 
single  quantum  well  (SQW)  GaAs  lasers  will  show  characteristic  temperature 
(Tq)  discontinuities  accompanied  by  large  jvunps  in  operating  wavelength 
[Zory  86].  Despite  the  complexity  of  these  phenomena,  it  has  been  possible 
to  model  them  successfully  by  assuming  that  the  optical  gain  is  derived  from 
k-conserved  stimulated  processes  between  free  electrons  and  holes  in  quan- 
tum well  subbands.  As  a consequence,  many  people  now  believe  that  this 
type  of  EHP  model  is  the  correct  one  to  use  when  designing  III-V  QW  lasers. 

However,  there  have  been  reports  that  the  optical  gain  in  CdZnSe  QW 

« 

lasers  is  not  derived  from  stimulated  processes  between  free  electrons  and 
holes  but  between  bound  electron  and  holes  or  excitons  [Ding  90,  Ding  92, 
Pele  92].  As  evidence,  measurements  of  absorption  spectra  of  quasi-three 

O 

dimensional  (active  layer  thickness  d=200A)  and  two  dimensional  (active 

o 

layer  thickness  d=30A)  samples  at  various  temperatures  were  used  [Ding  90]. 
It  has  been  observed  that  the  dominant  exciton  peak  at  low  temperature  dis- 
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o 

appeared  at  about  200K  for  d=200A  sample  while  the  n=l  HH  transition 
peak  remains  distinct  at  temperatures  well  beyond  room  temperature  for 
d=30A  sample.  However,  it  should  be  noted  that  the  step  like  nature  of  the 
density  of  state  of  a QW  structure  could  give  rise  to  such  an  absorption  peak 
which  can  be  observed  even  above  room  temperatures.  In  fact,  those  peaks 
have  also  been  observed  for  GaAs  QW  structures,  and  hence  it  is  questionable 
if  an  “excitonic  gain  (EX)”  model  rather  than  the  EHP  model  would  be  appro- 
priate for  designing  CdZnSe  QW  lasers.  As  shown  in  the  previous  section, 
the  EHP  model  can  be  used  to  make  quite  accurate  predictions  of  the  temper- 

m 

ature  dependence  of  for  CdZnSe  QW  lasers.  To  further  investigate  the 
ability  of  the  EHP  model  to  predict  QW  lasing  phenomena  in  CdZnSe  lasers, 
we  decided  to  see  if  we  could  obtain  dual  wavelength  operation  (simultaneous 
n=l  and  n=2  transition  lasing).  In  order  to  make  sure  that  our  predictive 
capabilities  and  experimental  techniques  were  correct,  we  decided  to  first 
obtain  dual  wavelength  operation  on  AlGaAs  QW  lasers. 

5.2  Dual  Peak  Wavelength  Operation  of  AlGaAs  QW  Lasers 

For  GaAs  QW  lasers,  it  has  been  shown  that  n=2  transition  lasing  can 
be  achieved  by  increasing  cavity  losses  siifficiently  [Zory  86].  There  are  a 
number  of  ways  to  increase  the  cavity  losses;  make  the  cavity  length  short, 
increase  the  device  temperature,  or  make  a stripe  width  narrower,  etc. 
Utilizing  short  cavity  lasers,  a thermally  induced  peak  wavelength  shift  of 
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about  500A  was  observed  from  50A  AlGaAs  QW  diode  lasers.  In  that 
experiment,  the  peak  wavelength  shift  occurred  between  58 °C  and  60°C.  In 
order  to  make  the  wavelength  switching  occur  at  low  carrier  densities  at 
room  temperature,  a thicker  AlGaAs  QW  laser  (thickness=160A)  was  chosen 
as  shown  in  Fig.  5.1.  The  energy  dispersion  was  calculated  using  multiband 
effective  mass  theory  (k  • p method)  for  the  valence  band  and  the  parabolic 
approximation  for  the  conduction  band.  The  results  are  shown  in  Fig.  5.2. 
In  order  to  increase  the  cavity  losses,  short  cavity  length  lasers  with  various 
stripe  widths  were  used.  The  spectral  gain  and  spontaneous  emission  rate 
was  calculated  using  the  EHP  model  discussed  previously.  The  resialts  are 
shown  in  Fig.  5.3  and  Fig.  5.4.  An  important  feature  can  be  found  in  the 
calculated  results;  at  some  injection  level,  N=3X18/cm^  for  instance,  gain  is 
mayirrmm  at  a low  photon  energy  (n=l  transition)  while  the  spontaneous 
emission  rate  peaks  at  a high  photon  energy  (n=2  transition).  It  means  that 
even  though  spontaneous  emission  at  n=2  is  dominant  at  this  carrier  density, 
lasing  action  will  occurr  at  n=l.  This  can  be  explained  by  appreciating  the 
formula  for  the  spectral  gain  and  the  spontaneous  emission  rate  given  in  Eq. 
(4.2)  and  Eq.  (4.4),  respectively.  In  terms  of  the  thermal  carrier  distribution, 
the  spontaneous  emission  rate  depends  on  the  multiplication  of  the  electron 
and  hole  probability  distribution  functions,  f^d-fy).  However,  the  spectral 
gain  depends  on  the  difference  between  the  distribution  functions,  fg-fy.  On 
the  other  hand,  at  N=4.4X18cm^,  the  peak  gain  for  n=l  and  n=2  are 
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Figure  5.1  Epitaxial  structure  of  AlGaAs  160 A QW  lasers 
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Figure  5,2  The  energy  vs.  in-plane  k vector  dispersion  for  160A 

AlGaAs  QW  lasers. 
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Figure  5.3  The  spectral  gain  of  AlGaAs  160A  QW  lasers  for  various 

injection  current  at  300K. 
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Figure  5.4  The  spontaneous  emission  rate  of  AlGaAs  160A  QW  lasers 

for  various  injection  current  at  300K. 
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comparable.  If  the  peak  gain  at  this  carrier  density  (about  1350cm'^)  is 
about  equal  to  the  threshold  gain,  then  it  might  be  possible  to  see 
simultaneous  lasing  at  n=l  and  n=2  wavelengths.  Hence,  the  diode  laser 
cavity  length  was  determined  so  that  g^jj  = 1350cm‘^  from  Eq.  (4.12).  For 
this  AlGaAs  QW  structure,  F =0.06,  Ri=R2=0.32,  aj=5.7cm'^.  So  L=150pm 
was  chosen  and  gth=1360cm'^.  The  measured  -spectra  for  stripe  width(w) 
10|im,  25|im,  and  50|im  at  various  carrier  densities  are  shown  in  Fig.  5.5,  Fig. 
5.6,  and  Fig.  5.7,  respectively.  As  can  be  seen  in  the  figures,  it  is  quite  clear 
that  both  the  n=l  and  n=2  peak  wavelengths  compete  at  low  current  injection 
level.  As  the  current  is  increased,  one  of  the  transitions  becomes  more 
distinct.  At  high  injection,  lasing  action  was  observed  at  775nm  for  narrow 
stripe  lasers  (w=10nm),  while  for  the  wide  stripe  lasers,  lasing  occurred  at 
SOOnm.  This  peak  lasing  wavelength  difference  (25nm)  corresponds  to  the 
energy  difference  between  the  maximum  peak  gains  for  n=l  and  n=2  shown 
in  Fig.  5.3.  The  peak  wavelength  difference  between  the  wide  and  the 
narrow  stripe  lasers  can  be  explained  in  terms  of  lateral  (along  the  junction 
plane)  mode  confinement.  In  the  estimation  of  the  confinement  factor  F, 
only  transverse  (perpendicialar  to  the  jimction)  confinement  was  considered. 
It  is  a good  assumption  that  the  field  intensity  is  well  confined  in  the  lateral 
direction  for  wide  stripe  lasers,  but  not  for  a narrow  stripe  laser  in  a gain 
guiding  structure.  For  narrow  stripe  lasers,  there  can  be  a sizable  mode  loss 
in  the  lateral  direction  and  hence,  effectively,  F in  Eq.  (4.12)  decreases.  As  a 
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Figure  5.5 
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The  measured  spectriim  of  AlGaAs  160A  QW  lasers  for 


L=150|im,  w=10pm  at  300K. 
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Figure  5.6 


The  measured  spectrum  of  AlGaAs  160A  QW  lasers  for 


L=150|im,  w=25  |im  at  300K. 
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Figure  5.7  The  measured  spectrum  of  AlGaAs  160A  QW  lasers  for 

L=150|im,  w=50|im  at  300K. 
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consequence,  for  a fixed  threshold  mode  gain  more  peak  gain  is  required 
to  compensate  the  reduced  effective  confinement  factor  for  narrow  stripe 
lasers.  The  cavity  length  was  decided  so  that  both  n=l  and  n=2  peak  gains 
reach  gjj^  simtdtaneously  for  wide  stripe  lasers  as  is  the  case  for  w=25nm  and 
w=50nm.  For  narrow  stripe  lasers,  larger  peak  gain  is  required  in  order  to 
compensate  the  lateral  mode  loss  and  such  gain  is  achieved  by  n=2  transition 
as  can  be  seen  in  Fig.  5.3. 

5.3  Dual  Peak  Wavelengths  of  CdZhSe  QW  Lasers 

In  order  to  reproduce  a similar  experiment  for  ZnSe  based  QW  lasers, 
thinner  QW  thickness  (100 A instead  of  160 A)  was  selected  because  the  effec- 
tive masses  of  ZnSe  are  heavier  than  those  of  GaAs,  and  hence  the  spacing 
between  each  quantized  energy  levels  is  smaller  for  same  thickness.  The 
effective  masses  of  GaAs  and  ZnSe  are  shown  in  Table  4.1.  The  laser  mate- 
rial used  in  this  experiment  was  grown  by  3M  Photonics  Research  Laboratory 
using  the  MBE  technique.  The  epitaxial  structure  is  identical  with  their 
first  demonstrated  blue-green  laser  as  shown  in  Fig.  5.8.  The  energy  disper- 
sion was  calculated  with  the  multiband  effective  mass  theory  and  the  results 
are  shown  in  Fig.  5.9.  The  spectral  gain  was  calculated  based  on  the  energy 

dispersion  for  80K,  150K,  and  300K  as  given  in  Fig.  5.10,  Fig.  5.11.  and  Fig. 

« 

5.12.  The  carrier  density  in  the  calculation  was  selected  so  that  both  n=l 
and  n=2  peak  gain  have  about  same  value.  The  experiment  was  restricted  to 
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Table  5. 1 

The  effective  masses  and  the  Luttinger  parameters  of  GaAs  and  ZnSe 


GaAs 

ZnSe 

nic 

0.067nio 

0.17nio 

<» 

mhh 

0.59mo 

0.78mo 

mih 

0.09mo 

0.145mo 

yi 

6.85 

3.77 

yi 

2.1 

1.24 

^3 

2.9 

1.67 
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p+  ZnSe 

p ZnSSe 

1.5pm 

p ZnSe 

1.0pm 

«8S8888ffl8888888888888iS888888888}B888fflffi^^ 

n ZnSe 

1.0pm 

n ZnSSe 

2.5pm 

n+ZnSe 

« 

n GaAs  Buffer 

n GaAs  substrate 

! M 1 iii 

Au  Electrode 
Pol3dmide 


Cd  oZn  oSe 
idoA  QW 


Indium  Electrode 


Figure  5.8  The  epitaxial  structure  of  ZnSe  lOOA  QW  lasers. 
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Figure  5.9  The  energy  dispersion  of  conduction  and  valence  band  of 

Cdo.2Zno,8Se  lOOA  QW  laser 
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Figure  5.10  The  convolved  and  unconvolved  gaiii  spectra  of  ZnSe 

lOOA  QW  lasers  for  N=4.5el8  cm‘^  at  80K. 
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Figure  5.11  The  convolved  and  unconvolved  gain  spectra  of  ZnSe 

lOOA  QW  lasers  for  N=6.5el8  cm"^  at  150K. 
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O 

Figure  5.12  The  convolved  and  imconvolved  gain  spectra  of  ZnSe  lOOA 

QW  lasers  for  N=9.8el8  cm'^  at  300K. 
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low  temperature  for  the  following  reasons;  First,  the  liaminescence  efficiency 
of  the  material  is  poor  and  it  did  not  lase  at  room  temperature.  Secondly,  the 
carrier  scattering  time  is  much  shorter  at  room  temperature  as  discussed  pre- 
viously in  section  4.3.  The  resulting  wider  line  shape  function  tends  to 
smooth  the  spectral  gain  curve  to  the  point  where  the  n=l  and  n=2  transi- 
tions cannot  be  resolved.  In  a similar  procedure  with  AlGaAs  160A  QW 
lasers  discussed  previously,  CdZnSe  diode  lasers  were  designed  so  that  both 
the  n=l  and  n=2  gain  peaks  have  the  same  value  at  threshold.  From  Fig. 
5.10  and  Fig.  5.11,  the  proper  gj^,  was  foimd;  2200cm‘^  for  150K,  and  3000cm' 
^ for  80K.  These  threshold  gains  at  each  temperature  can  be  achieved  by 
designing  the  cavity  length  accordingly.  The  corresponding  cavity  lengths 
are  430pm  for  80K  and  650pm  for  150K.  The  parameters  used  in  the  estima- 
tion are  a-=8.1cm'^,  Ri=R2=0.23  and  F =0.014.  Laser  diodes  were  fabricated 
using  the  following  procedure: 

Cleaning 

1. TCE  soft  boiling  2 minutes 

2. Acetone  rinse  2 times 

3.  Methanol  rinse  2 times 

4. Nitrogen  blow 
Photolithography  I 

l.P.R  Coating  5500  rpm  30  second  using  Photoresist  AZ1400-33 


2.Soft  Bake  90C  30  minutes 
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3.  Mask  align  and  expose  at  20mW/cm^  for  25  seconds 

4. Develop  for  30  seconds  with  developer  AZ  319 
Metal  Deposition 

o 

Deposit  lOOOA  of  gold  using  electron  beam  evaporator 
Metal  Lift-Off 

Acetone  Rinse 

Methanol  Rinse 
Photolithography  II 

l.P.R  Coating  5500  rpm  30  second  with  Photoresist  AZ 1400- 17 

2.Soft  Bake  at  90C  for  30  minutes 

3. Mask  align  and  expose  at  20mW/cm^  for  15  seconds 

4. Develop  for  20  seconds  in  the  developer  AZ  319 
Metal  Deposition 

Deposit  Ni  lOOOA/Au  lOOOA  using  electron  beam  evaporator 
Cleaving 

Cleave  the  wafer  with  250|im  cavity  length 

The  diode  laser  chips  were  scribed  from  the  laser  bars  and  attached  on 
copper  blocks  and  wire  bonded.  The  near  field  pattern  was  measured  for 
each  sample  and  poor  devices  screened  out.  The  samples  were  placed  in  a 
dewar  chamber  and  the  air  pressure  reduced  with  a cryogenic  prnnp  down  to 
30millitorr.  Then  liquid  nitrogen  was  poured  into  the  dewar  and  the 
chamber  connected  with  a temperature  controller  imit.  The  spectral 
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measurements  which  showed  the  most  evidence  for  dual  wavelength 
operation  are  shown  in  Fig.  5.13  and  Fig.  5.14.  Although  a dip  can  be  seen  in 
each  spectra,  the  data  was  recorded  at  sub-lasing  threshold  currents.  At 
threshold  and  above,  dual  wavelength  lasing  of  the  t}rpe  demonstrated  for  the 
AlGaAs  QW  lasers  could  not  be  observed.  The  reason  is  due,  we  believe,  to 
the  fact  that  carrier  scattering  process  in  ZnSe  based  materials  are  much 
faster  than  in  GaAs  based  materials. 

5.4  Summary 

Spectral  measurements  were  made  for  AlGaAs  and  CdZnSe  QW  lasers. 
The  QW  thickness  and  cavity  length  were  chosen  so  that  both  n=l  and  n=2 

Q 

peak  gain  reached  threshold  at  the  same  carrier  density.  For  AlGaAs  160A 
QW  lasers,  dual  wavelength  lasing  operation  was  clearly  seen.  For  CdZnSe 

o 

lOOA  QW  lasers,  those  peaks  were  not  quite  distinct  due  to  the  shorter  scat- 
tering time. 


Intensity  (a.u) 
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Figure  5.13  The  measured  spectrum  of  ZnSe  100 A QW  lasers  for  vari- 
ous injection  current  at  80K.  (L=430pm,  w=10|im) 


Intensity  (a.u) 
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Wavelength  (nm) 


Figure  5.14 


o 

The  measured  spectnun  of  ZnSe  lOOA  QW  lasers  for  vari- 
ous injection  current  at  150K.  (L=650|im,  w=25|im) 


APPENDIX 


For  k=0, 


H|n,  0)  = E„o|n,  0) 

where  EnO  is  the  eigenvalue  of  band  n at  k=0.  From  Eq.  (2.3), 

|n,  k>  = e**'  *'|n,  0> 

Hence, 
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Using  H from  Eq.  (2.1)  we  have 
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Hence, 
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where  is  defined  in  Eq.  (2.5). 
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